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I. INTRODUCTION 


Modern cosmology is founded upon the so-called 
Copernican principle, which posits that the universe 
‘looks’ on average to be the same regardless of where 
one is in the universe or in which direction one looks. In 
other words, every point in the universe is identical; no 
point is special. More formally, the Copernican princi¬ 
ple states that Cauchy surfaces of the universe can be 
admitted that are homogeneous and isotropic, and this 
symmetry can be expressed mathematically by writing 
the universe’s metric in the form 
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where aft ) is a time-dependent function known as the 
scale factor and k is a curvature constant. The sign of 
k determines whether Cauchy surfaces of constant time 
t. will be open, flat, or closed, with k < 0 being open, 
k = 0 being flat, and k > 0 being closed. The metric © 
is known as the Friedmann-Lemaitre-Roberts on- Walker 
(FLRW) metric. Though the form of this FLRW metric is 
fixed by consideration of symmetries alone, the function 
for the scale factor aft) is instead determined by general 
relativity. Inserting the metric into the Einstein field 
equations yields 



* R.Liu@damtp.cam.ac.uk 
1 R.M.Williams@damtp.cam.ac.uk 


where p and p are the energy density and pressure of 
any perfect fluid filling the space, and where A is the 
cosmological constant. This pair of differential equations 
is known as the Friedmann equations, and their solution 
determines a(t). 

These FLRW models have had great success in explain¬ 
ing much of the universe’s behaviour, including most no¬ 
tably the Hubble expansion of the universe, the cosmic 
microwave background, and baryon acoustic oscillations. 
Indeed, the underlying assumption of homogeneity and 
isotropy appears well-supported by precision measure¬ 
ments showing the CMB to be isotropic to within one 
part in 100,000 pj. Yet in spite of this, observations also 
clearly show that the late, matter-dominated universe is 
not homogeneous and isotropic except at the coarsest of 
scales. Instead, matter is distributed predominantly in 
clusters and superclusters of galaxies with large voids in 
between, and the physical effects of such a ‘lumpy’ uni¬ 
verse are still not fully understood. 

Indeed, there has been intense interest recently over 
the possible importance of inhomogeneities to observa¬ 
tional cosmology. Perhaps the area of greatest interest 
concerns the possible effects of inhomogeneities on recent 
redshift measurements from Type la supernovae (SNla). 
When fitted to perfectly homogeneous FLRW models, 
these measurements have led to the conclusion that the 
universe’s expansion is accelerating [g], and to account 
for this acceleration, cosmologists have posited the ex¬ 
istence of some exotic matter, known generally as dark 
energy. However other cosmologists have instead posited 
that much, if not all, of this acceleration is as an apparent 
effect, arising from fitting data from an inhomogeneous 
universe onto a homogeneous model 0. It has been ar¬ 
gued that this effect actually arises from the non-linear 
structure of the inhomogeneities and that such a struc¬ 
ture could not be adequately modelled simply by per- 
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turbing an FLRW universe 0- Therefore, as there has 
so far been no confirmed direct observation of any exotic 
matter, the relative importance of dark energy compared 
to inhomogeneities in explaining the supernovae data re¬ 
mains an open question. 

For this reason, there has been a resurgence of inter¬ 
est recently in the various approaches to building non- 
perturbative, inhomogeneous models. One notable ex¬ 
ample is the Lemaitre-Tolman-Bondi (LTB) models. It 
has been shown that such models can account for all 
known cosmological observations, including SNla red- 
shifts, without requiring any exotic matter, but the ob¬ 
server must sit at the centre of a Hubble-scale under- 
dense region Q. In other words, the observer sits at the 
centre of an isotropic but non-homogeneous universe. 

On the other hand, the ‘Swiss-cheese’ models of Ein¬ 
stein and Straus retain an FLRW background but intro¬ 
duce inhomogeneities by replacing co-moving spherical 
FLRW regions with Schwarzschild or LTB regions @, 7|. 
By an appropriate fitting of these regions into the FLRW 
background, the resulting space-time will still be an ex¬ 
act solution to the Einstein field equations. However, 
because of their FLRW background, these models will 
still be dynamically identical to FLRW universes. And 
although there is no a priori reason to believe the optical 
properties should be identical as well, recent studies have 
shown that they are, in fact, broadly similar E}. 

In this paper, we shall consider a different class of uni¬ 
verses, one that is not in any way based on FLRW uni¬ 
verses. We shall consider the so-called lattice universes 
where the matter content on each Cauchy surface consists 
of identical point masses arranged into a regular lattice. 
Unlike the LTB or Swiss-cheese models, the lattice uni¬ 
verse has a truly discrete matter content and is other¬ 
wise vacuum throughout, which is more representative 
of the actual universe’s matter distribution. Because of 
the regular arrangement of the point masses, these lattice 
universes still possess a high degree of symmetry, though 
not as great as that of FLRW universes. 

We shall focus on lattices formed by tessellating 3- 
spaces of constant curvature with identical regular poly¬ 
hedral cells. The possible lattices that can be constructed 
from such a tessellation have been summarised in Ap¬ 
pendix [A] and we shall refer to these lattices as Coxeter 
lattices. To ‘construct’ a lattice universe then, we select 
one of the Coxeter lattices of Appendix lAl and distribute 
a set of point masses in a regular manner on that lat¬ 
tice, such as at the centres of the cells, at the centres 
of the faces, at the mid-points of the edges, or at the 
vertices. Naturally, after the masses have been arranged 
in this manner, the Cauchy surfaces will no longer be 
surfaces of constant curvature; however, the metric is 
still expected to be invariant under the same symmetry 
transformations that leave the lattice invariant, symme¬ 
tries which include discrete translation symmetries, dis¬ 
crete rotational symmetries, and reflection symmetries at 
the cell boundaries. In other words, the lattice universe 


should have a metric of the form 

ds 2 = — dt 2 + 7 ^ ( t , x ) dx a dx b , (4) 

where (t , x) is the 3-dimensional metric for constant 

t hypersurfaces, and Latin indices a, b = 1, 2, or 3 denote 
spatial co-ordinates only; the spatial metric 7 ® (f, x ) at 
constant t would possess the lattice symmetries. Effec¬ 
tively, the Copernican symmetries of FLRW universes 
have been reduced to just these symmetries. This paper 
will focus exclusively on closed lattice universes based on 
the tetrahedral Coxeter lattices; these lattices consist of 
5, 16, or 600 identical, equilateral tetrahedral cells. 

There has been a variety of approaches adopted to 
modelling such universes. Using exact methods, Wheeler 
[Sj as well as Clifton et al. p] have successfully con¬ 
structed the exact 3-metric (t = 0,®) for the time- 

symmetric Cauchy surfaces of closed universes. Ko- 
rzynski has further generalised this work by examining 
the case where there is an arbitrary number of masses, 
not necessarily arranged in a lattice M- Clifton et 
al. have also examined the dynamics of the closed 
universes by evolving their initial data along certain 
highly-symmetric curves Bruneton and Larena have 
modelled the dynamics of the flat universe using an 
exact but perturbative expansion of the metric about 
Minkowski space-time fl2| . Using numerical approaches, 
Bentivegna and Korzynski have studied the evolution of 
closed universes from initial data on a hypersurface at 
time-symmetry (T3j . while Yoo et al. as well as Ben¬ 
tivegna and Korzynski have studied the dynamics of the 
flat universe EH EH- On the other hand, Lindquist and 
Wheeler [u| have devised an approximation to the lat¬ 
tice universe, generalised by Clifton and Ferreira E3, 
wherein each polyhedral lattice cell gets approximated by 
a spherical cell with Schwarzschild geometry inside. This 
has been applied to study the evolution of the closed, 
flat, as well as open universes EMI. 

In this paper, we shall consider another approach to 
modelling the lattice universe; we shall adopt a Regge 
calculus formalism originally developed by Collins and 
Williams (CW) to model closed FLRW universes El- 
Regge calculus [2C| is a highly versatile formalism that 
can in principle approximate any solution of the Ein¬ 
stein held equations using a piece-wise linear manifold; 
this makes it particularly suitable to studying systems 
where an exact solution is difficult to obtain. The Regge 
manifold is constructed by gluing hat blocks together 
such that neighbouring blocks share an entire face; as 
the blocks are hat, the metric inside is the Minkowski 
metric. The Regge manifold is generally referred to as a 
skeleton. The solutions of Regge calculus are generally 
expected to converge at second order in the skeletal edge- 
lengths to the corresponding continuum solutions of the 
Einstein held equations [2l| ■ Although we shall focus on 
closed lattice universes in this paper, much of the work 
is readily generalisable to hat and open universes. 

Collins and Williams have constructed their skeleton 
from a one-parameter family of space-like Cauchy sur- 
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faces that foliate the entire skeleton. Each Cauchy sur¬ 
face triangulates a closed FLRW Cauchy surface with 
equilateral tetrahedra such that all vertices, edges, and 
faces in the triangulation are identical. We note that 
closed FLRW Cauchy surfaces of constant t can be em¬ 
bedded as 3-spheres in 4-dimensional Euclidean space 
E 4 . The scale factor a(t) can always be re-scaled so that 
the FLRW curvature constant k becomes unity, in which 
case, a(t) would equal the 3-sphere radius of the embed¬ 
ding. The embedding is then given by 

r = sin x, 
x 1 = a(t) cos x, 

x 2 = a(t) sin x cos 0, (5) 

x 3 = a(t ) sin x sin 9 cos </>, 
x 4 = a(t) sin x sin 9 sin (f>, 

for 0 < Xi 0 < tt and 0 < 4> < 27r. CW Cauchy surfaces 
would triangulate such 3-sphere Cauchy surfaces with 
identical equilateral tetrahedra, and according to Coxeter 
HH, such a triangulation is only possible using 5, 16, and 
600 tetrahedra. In fact, the triangulations correspond to 
the three closed tetrahedral Coxeter lattices of Appendix 
m this makes the CW skeletons particularly apposite for 
modelling lattice universes, as their Cauchy surfaces nat¬ 
urally provide suitable lattices into which the masses of 
the lattice universe could be embedded. Table Q] tabulates 
the numbers of vertices, edges, triangles, and tetrahedra 
for each of the possible triangulations. The CW Cauchy 

TABLE I. The number of simplices in each of the three trian¬ 
gulations of the 3-sphere with equilateral tetrahedra as well 
as the number of triangles meeting at any edge. We introduce 
N 3 , N 2 , N\, and No to denote the number of tetrahedra, tri¬ 
angles, edges, and vertices in the Cauchy surface. 


Tetrahedra 

(Ns) 

Triangles 

(A 2 ) 

Edges 

(Ni) 

Vertices 

(No) 

Triangles 
per edge 

5 

10 

10 

5 

3 

16 

32 

24 

8 

4 

600 

1200 

720 

120 
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surfaces are then joined together by a series of time-like 
edges called struts connecting each vertex on one surface 
with its time-evolved image on the next. Because all ver¬ 
tices on a surface are identical, all struts between any two 
surfaces are identical as well. With this construction, the 
CW Cauchy surface at discrete time parameter ti can be 
characterised by just two distinct lengths, the tetrahedral 
edge-length l(ti) and the strut-length m(ti). Surfaces at 
different ti are completely identical apart from an overall 
re-scaling of the l(ti) length-scale. Therefore, Collins and 
Williams interpreted Z (i*) to be the Regge calculus ana¬ 
logue of a(t). By then solving the Regge action and then 
taking the continuum time limit (t i+ 1 — ti) —> 0, they 
were able to derive an expression for the length-scale as 
a function l(t) of continuum time t. The CW skeletons 


were first applied to model closed dust-filled FLRW uni¬ 
verses and were found to yield very accurate results, with 
the models with agreater number of tetrahedra yielding 
higher accuracy [ 13 , [23 • 

In this paper, we shall also consider a different type 
of lattice universe. We shall consider a closed lattice 
universe where a single mass gets perturbed in magni¬ 
tude. However as a prelude to this investigation, we 
shall first explore the properties of the completely un¬ 
perturbed lattice universe. Thus, this paper is organised 
as follows. We shall begin, in the second section, with 
a brief exposition of Regge calculus with CW skeletons. 
In the third section, we then explore the behaviour of 
the unperturbed lattice universe. We shall show that the 
universe’s behaviour depends on where the masses are lo¬ 
cated: the universe becomes unconditionally divergent if 
the masses are placed at the vertices of the Cauchy sur¬ 
face but is unconditionally convergent if the masses are 
placed anywhere in a spherical region around the centres 
of the tetrahedra. Thereafter, we shall work only with 
models where the masses are at the tetrahedral centres. 
We also make this choice for the following reason: in 
FLRW universes, a test particle that is co-moving with 
the universe would be following a geodesic as well, and 
we suspect this to also be the case in the lattice universe; 
thus if we want particles in our model also to be both co¬ 
moving and following geodesics across the entire Regge 
space-time, then we must place the particles at the cen¬ 
tres of the tetrahedra; at any other location, co-moving 
particles will not follow global geodesics. In the fourth 
section, we shall perturb one of these masses, construct 
the corresponding perturbed Regge model, and derive the 
relevant equations governing the model’s evolution. We 
shall focus exclusively on the 5-tetrahedra model: this 
model would involve only two sets of perturbed edges 
whereas the other two models would require many more, 
with each set having its own independent length. In the 
fifth section, we shall consider the application of the ini¬ 
tial value equation at the moment of time symmetry to 
this model: we shall derive certain conditions that the 
initial conditions of the Regge equations must satisfy in 
order to be consistent with this equation. In the final 
section, we shall examine the behaviour of the model for 
various perturbations and compare it against that of the 
unperturbed model; we shall then close with a brief dis¬ 
cussion of certain assumptions inherent in our model. 

In this paper, we shall use geometric units where G = 
c = 1. 


II. REGGE CALCULUS WITH CW 
SKELETONS 

As mentioned above, Regge calculus approximates any 
continuous space-time using a piece-wise linear manifold 
composed of flat blocks glued together at their faces. 
Regge calculus customarily uses 4-simplices as the blocks, 
though as we shall soon see with the CW skeleton, this 
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is not always the case. 

Curvature in the skeleton manifests itself as conical sin¬ 
gularities concentrated on the sub-faces of co-dimension 
2; these sub-faces are known as hinges. Each block will 
have two sub-faces of co-dimension 1 meeting at any of 
its hinges; we shall refer to such sub-faces simply as the 
faces of the block. If a hinge were flat, then the dihedral 
angles between all faces meeting at the hinge would sum 
to 27 t; any deviation from 2-7 t provides a measure of the 
curvature and is known as the deficit angle. The deficit 
angle <5,; at a hinge labelled i is given by 

k = 2ti-- ( 1 * * * * 6 ) 

0 

where 6j is the dihedral angle at the hinge formed by 
the faces of the block labelled j and the summation is 
over all blocks meeting at the hinge. 

If the skeleton consists solely of 4-simplices, then to 
completely specify the skeletal geometry, one need only 
specify the lengths of all edges. This follows because the 
internal geometry of any n-simplex, including its angles 
and the areas of its sub-simplices, is completely deter¬ 
mined when the lengths of its C(n+ 1, 2) edges are spec¬ 
ified. Thus the skeletal edge-lengths serve as the Regge 
analogue of the metric; in analogy to how the metric 
is determined by the Einstein field equations in general 
relativity, the edge-lengths are determined by the Regge 
field equations, the Regge analogue to the Einstein field 
equations. 

The Regge field equations are obtained by following 
a variational approach similar to how the Einstein field 
equations are obtained from the Einstein-Hilbert action. 
In general relativity, the Einstein field equations can be 
derived by varying the Einstein-Hilbert action 

S EH = ^j R V=9 di x - Y, M 'J ds i ( 7 ) 

i £ {particles} 

with respect to the metric tensor g where R is the Ricci 
scalar, g = det(g^), Mi is the mass of particle i, and dsi 
is its line element; the summation is over all particles in 
the space-time. When applied to a Regge skeleton, this 
reduces to the Regge action [2Cj 

SRegge = Ai Si — Mj Sjj, (8) 

i £ {hinges} i £ {particles} 
j £ {blocks} 

where Ai is the area of a hinge in the Regge skeleton, 
Si its corresponding deficit angle, and the length of 
particle i’s path through block j; the first summation is 
over all hinges in the skeleton while the second is over 
all particles and all blocks of the skeleton. Note that 
if particle i never passes through block j, then Sij will 
accordingly be zero. 

Since the skeletal edge-lengths are the Regge analogue 
of the metric, the Regge action is varied with respect to 


an edge-length £j to get the Regge field equations 

0= w 

z£ {hinges} J i £ {particles} J 
j £ {blocks} 

where the variation of the deficit angles has cancelled out 
owing to the well-known Schlafli identity [2C| , 


this identity holds for any individual block k, with the 
summation being over all hinges in the block and 6^ 
being the block’s dihedral angle at hinge 

Collins and Williams originally constructed their skele¬ 
ton not out of 4-simplices but out of 4-blocks; each 4- 
block corresponds to the truncated world-tube of a tetra¬ 
hedron as it evolves from one Cauchy surface to the next. 
If we denote the Cauchy surface at time ti by £* and the 
surface at time U+\ by Ej+i, then a 4-block between £* 
and £ l+ i consists of a tetrahedron in £j with edges of 
length U = l(ti ), a tetrahedron in £ i+ i with edges of 
length li -|_i = l(ti- |_i), and four equal-length struts con¬ 
necting each tetrahedral vertex in to its time-evolved 
counterpart in S^+i. Figure Q] depicts a typical 4-block; 
for simplicity, we shall sometimes refer to the tetrahedron 


O' 



FIG. 1. An equilateral tetrahedron of edge-length U at time 
ti evolves to a tetrahedron of edge-length h+i at time ti+i, 
tracing out a 4-dimensional world-tube. The struts are all of 
equal length. 


1 In the standard formulation of Regge calculus, one actually uses 

a simplicial manifold where every block is a 4-simplex, and the 
Schlafli identity is usually formulated in terms of simplices rather 

than arbitrary blocks. However, any block can always be trian¬ 

gulated into simplices, and one can then apply the simplicial 

form of the Schlafli identity to the triangulated block to obtain 

the form of the identity we have above, using the chain rule if 

necessary to satisfy any constraints on the block’s geometry. 
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in Ej+i as the upper tetrahedron and the tetrahedron in 
Si as the lower tetrahedron, as that is how they appear 
in the figure. 

Since we are no longer using 4-simplices as the skele¬ 
tal building block, specification of the edge-lengths will 
in general not be sufficient to determine the entire skele¬ 
tal geometry. There are two approaches one could take 
to completing the geometry. The first is to fully trian¬ 
gulate each 4-block into 4-simplices by introducing ex¬ 
tra edges; one would then have to specify the lengths 
of these new edges. The second, which was taken by 
Collins and Williams, is to specify the 4-block’s inter¬ 
nal geometry. Regardless of which approach is taken, 
the 4-block’s geometry would be determined by two re¬ 
quirements: (i) that all struts have the same length; 
and (ii) that there be no twist or shear along the 4-block. 
These requirements would manifest themselves either as a 
specification of the new edge-lengths in the first approach 
or as a specification of the 4-block’s internal geometry in 
the second. Brewin has likened the requirements to a 
choice of lapse and shift function in the ADM formal¬ 
ism. Indeed, the standard form of the FLRW metric © 
also implies a certain foliation of FLRW space-time, and 
Collins and Williams’ choice seems closest to the lapse 
and shift implicit in this foliation. 

These constraints imply that the lower tetrahedron 
would simply expand or contract uniformly about its cen¬ 
tre when evolving to the upper tetrahedron. The geome¬ 
try can best be understood by introducing a co-ordinate 
system for the 4-block. As depicted in Figure [I] we label 
the vertices of the lower tetrahedron by A, B , C, D, and 
their counterparts in the upper tetrahedron by A', B ', 
C', D ', respectively. The lower vertices’ co-ordinates are 
then 


A 

B 

C 

D 


li li 


U 


2 ’ 2 ^ 3 ’ 2V% 

li li li 
2’ 2 v / 3’ 2-\/6 


it, 


it. 


o,4,- 


k 


V 3 ’ 2%/6 


, iti \, 




( 10 ) 


the co-ordinates of their upper counterparts are given 
by an analogous expression where each letter becomes 
primed and each subscript changes from * to i + 1. Al¬ 
though a Euclidean metric is being used, the imaginary 
unit l has been introduced to the time co-ordinate so 
that inner products would effectively yield a signature 
of (+,+,+,—). In Collins and Williams’ approach, the 
4-block’s internal geometry is constrained to be that rep¬ 
resented by these co-ordinates for given edge-lengths k, 
and uii , that is, where the tetrahedron simply ex¬ 
pands or contracts uniformly about its centre. 

If we choose to fully triangulate the skeleton instead, 
then we shall introduce diagonals of type AD', BD', 


CD', AC', BC ', and AB' in each 4-block; this corre¬ 
sponds to one diagonal above each of the tetrahedral 
edges, as illustrated in Figure O and these diagonals di¬ 
vide the 4-block into four distinct 4-simplices, ABCDD', 
ABCC'D', ABB'C'D', and AA'B'C'D'. Each 4-block 



FIG. 2. The world-sheets generated by the six tetrahedral 
edges and their triangulation into triangular time-like hinges. 

in the skeleton will get triangulated in this manner, and 
it can be shown that this will lead to a consistent trian¬ 
gulation of the entire skeleton. To be consistent with the 
geometry described by m, the diagonals would then be 
constrained to have a length given by 

d? = -^lf + ^ (3 h+ i + k) 2 — Stf , (11) 

where Sti denotes the quantity Sti := tj + 1 — tj. 

When it comes to varying the Regge action, the stan¬ 
dard approach would be to regard each edge as being 
independent of all others; the action would then be var¬ 
ied one edge-length at a time, keeping all other lengths 
constant, and each edge would generally lead to an inde¬ 
pendent Regge equation. We shall refer to this approach 
as local variation. Local variation is always performed 
on a fully triangulated skeleton so as to avoid ambi¬ 
guities in the varied skeleton’s geometry. Consider for 
instance varying the edge AB in quadrilateral ABA'B' 
from length lf B to lf B : without any constraints on the 
quadrilateral’s internal geometry, there would be a wide 
range of possibilities that the varied quadrilateral could 
be, all possibilities having struts of length rrn and tetra¬ 
hedral edges of lengths lf B and lC[ . Thus constraints 
on the varied geometry are needed to determine a unique 
possibility. Under local variation, the constraint is im¬ 
posed by specifying a diagonal’s length and then keeping 
that length constant when another edge is varied. 

Collins and Williams however use a different approach 
to vary the skeleton: they would vary entire sets of identi¬ 
cal edges at once; for instance, they would vary all tetra¬ 
hedral edges on a Cauchy surface at once, or all struts 
between a pair of consecutive surfaces. Under this varia¬ 
tion, the tetrahedra in each 4-block would remain equilat¬ 
eral, the struts connecting them would continue having 
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equal length, and the 4-block itself would continue hav¬ 
ing no twist or shear; thus even when varied, all 4-blocks 
would continue having the same internal geometry as 
that specified by the 4-block co-ordinate system (flOl) . We 
shall refer to this approach as global variation. In stan¬ 
dard general relativity, global variation would be analo¬ 
gous to requiring the metric in the Einstein-Hilbert ac¬ 
tion to be of FLRW form (JTJ , and then varying the action 
with respect to a(t); this effectively imposes Copernican 
symmetries on the metric prior to varying the Einstein- 
Hilbert action. On the other hand, local variation is anal¬ 
ogous to the more standard approach where the Einstein- 
Hilbert action is varied first, yielding the Einstein field 
equations, and then the metric is set to be of FLRW form. 

However we have previously modelled the A-FLRW 
universe using the CW formalism [2ll | , and in those mod¬ 
els, we showed that locally varying the Regge action for 
a CW skeleton led to a set of non-physical equations. 
When we triangulated the skeleton, we had assumed that 
all tetrahedral edges in a Cauchy surface and all struts 
between a pair of surfaces would remain identical. Thus, 
when we locally varied the skeleton, we would set the 
relevant edges to be equal to help simplify the result¬ 
ing Regge equations; this has the effect of reducing the 
number of independent Regge equations to just one per 
distinct set of edges. We were effectively assuming that 
had we not set the lengths equal, the complete set of 
Regge equations would reveal them to be equal anyway. 
However, we believe this assumption to be unfounded be¬ 
cause we believe the diagonals actually disrupt the sym¬ 
metry between edges, rendering them no longer identical 
to each other. The diagonals are not actually distributed 
uniformly across the skeleton: some vertices are attached 
to more diagonals than others; for instance in our 4-block 
triangulation above, vertex A gets attached to three di¬ 
agonals while vertex D gets attached to none. Thus, 
the struts, which correspond to the vertices’ world-lines, 
would no longer be identical. For similar reasons, other 
geometric objects, such as tetrahedral edge-lengths and 
their associated world-sheets, may no longer be identical 
either. Therefore, it no longer makes sense to set edge- 
lengths to be equal after locally varying the skeleton; 
rather, each edge would have to be evolved independently 
of all others using its own Regge evolution equation, and 
there may potentially be one independent equation for 
each edge in the Cauchy surface. We did not encounter 
any similar problem when we globally varied the Regge 
action; in that case, all edges of the same type would in¬ 
deed be truly identical to each other. Thus drawing on 
this lesson, in this paper, we shall consider only solutions 
to the global Regge equation for all models. 

However, Brewin [231 ] has shown that, in certain cases, 
the global Regge equation can be related to the local 
one through a chain rule. In such cases, the local solu¬ 
tions would form a subset of the global ones. Although 
we have chosen to consider only global solutions, for the 
perturbed lattice universes, we shall derive the global so¬ 
lutions by using such a chain rule. This chain rule re¬ 


lationship between the global and local Regge equations 
will be further elaborated on when we consider the per¬ 
turbed universes in the fourth section. 

Brewin has also drawn several analogies between the 
ADM formalism and the CW formalism. He has likened 
the tetrahedral edge-lengths to the 3-metric of an ADM 
foliation and the Regge equations obtained from varying 
the tetrahedral edges to the ADM evolution equations. 
He has also likened the struts and diagonals to the ADM 
lapse and shift functions, respectively, and the Regge 
equations obtained from their variation to the ADM 
Hamiltonian and momentum constraints, respectively. 
Thus in this paper, we shall refer to the Regge equations 
obtained from the tetrahedral edges as evolution equa¬ 
tions, from the struts as Hamiltonian constraints, and 
from the diagonals as momentum constraints. In our 
study of the A-FLRW universe, we also found the Hamil¬ 
tonian constraint to be a first integral of the evolution 
equation; thus we could study the universe’s behaviour 
from the Hamiltonian constraint alone. Similarly, it can 
be shown that for the models of the unperturbed lattice 
universe, which will be considered below, the Hamilto¬ 
nian constraints are also first integrals of the evolution 
equations, though we shall not provide the proof. We 
shall assume a similar conclusion holds for models of the 
perturbed lattice universe. This conclusion means that 
we can determine the models’ evolution from their con¬ 
straint equations alone, and this is what we shall do. 

Before leaving this section, we wish to make a fi¬ 
nal comment on using co-ordinate system ms to cal¬ 
culate geometric quantities in the CW skeleton. This 
co-ordinate system greatly facilitates the calculation of 
any such quantity, but these quantities would get ex¬ 
pressed in terms of the time difference Sti rather than 
purely in terms of skeletal edge-lengths. In Regge calcu¬ 
lus, it is the skeletal edge-lengths that are to be varied, so 
we shall have to convert Sti into edge-lengths. Only the 
strut-lengths mi depend on Sti , since increasing the time 
separation between a pair of consecutive Cauchy surfaces 
lengthens to, but leaves U and U +1 unchanged. Therefore 
by using (in to calculate a strut-length, such as that of 
AA', we obtain the relation 

= (P’ 2 ~ ^ St >’ ( 12 ) 

where we have introduced the notation 



III. REGGE CALCULUS OF CLOSED, 
REGULAR LATTICE UNIVERSES 

We shall now apply the CW formalism to model 
closed, regular lattice universes. In this context, the 
Regge action ([8|) can be greatly simplified. First, a non- 
triangulated CW skeleton consists of only two distinct 
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types of hinges, space-like triangular hinges correspond¬ 
ing to the triangles of the equilateral tetrahedra and time- 
like trapezoidal hinges generated by the world-sheets of 
the tetrahedral edges as they evolve from one Cauchy sur¬ 
face to the next. An example of a triangular hinge would 
be ABC in the 4-block described by (ITUl) . and an exam¬ 
ple of a trapezoidal hinge would be ABA' B'. Secondly, 
in the regular lattice universe, all masses are identical. 
Hence, the action © can be expressed as 


87r5 = ^Af ap ^ rap + E^ tri -87riV p M £ *, 

. f trapezoidal 1 . f triangular 1 

hinges J hinges / 

(13) 

where M is the common mass of each particle, N p the to¬ 
tal number of particles in the universe, and s* the length 
of one particle’s trajectory between Cauchy surfaces £,; 
and £ i+ i. As we shall be varying with respect to the 
struts alone, the space-like triangular hinges can be ig¬ 
nored. 

When this action is varied globally with respect to the 
struts rrij. we obtain the Regge equation 


i 


a/(trap 

c trap 

drrij 1 


8nN p M E 

iefti} 


dsj 
drrij ’ 


(14) 


the general location of 

Vi = aA + /3B + 7C + iD, 

with vectors A, B, C, D denoting the position vectors 
of the tetrahedron’s four vertices relative to the tetra¬ 
hedral centre, and with constants a, (3, 7 , <5 satisfying 
0 < 01,13,7,5 < 1 and a + (3 + 7 + 5 = 1. Each mass 
is located at a distance |iq| from the tetrahedral centre 
given by 


Vi ■ Vi = i [3(a 2 + /3 2 + 7 2 + 5 2 ) 

O 

— 2 (a/3 + cry + a5 + /3y + f36 + 7 J)] l 2 , 

(15) 

with k being the length of a tetrahedral edge; as \v.j\/li 
is a constant of time, we shall denote this constant by v. 

We shall work with 4-block co-ordinates given by mni) 
and their £; + i counterparts to determine all geometric 
quantities; we can consider position vectors A, B, C, D 
as being identical to the position of vertices A , B , C, D 
in m■ In these co-ordinates, the length s, of each mass’ 
line element between £j and Ej+i is given by 

s 2 = v 2 5l 2 - St 2 , 


where the first summation is still over all trapezoidal 
hinges. 

There are several possible ways to arrange the masses 
into a regular lattice on a Cauchy surface; examples in¬ 
clude placing masses at the centres of the tetrahedra, 
the centres of the triangles, the mid-points of the tetra¬ 
hedral edges, or the tetrahedral vertices. Each of these 
configurations will yield a regular lattice; the new cell 
boundaries would lie along planes equidistant to pairs of 
masses that are nearest neighbours to each other, and the 
masses would consequently lie at the centres of the new 
cells. A 2-dimensional analogue has been illustrated in 
Figurc[9]of Appendix [Al where the original lattice, drawn 
in solid lines, consists of equilateral triangles tessellating 
flat 2 -dimensional space; masses placed at the mid-points 
of the triangular edges result in a new lattice consisting 
of rhomboidal cells, drawn in dashed lines. However as 
discussed in Appendix[A] not all new lattices would corre¬ 
spond to Coxeter lattices; those that do not would have 
non-regular polytopes as lattice cells and would conse¬ 
quently have reduced lattice symmetries. However we 
see that in this way, the CW formalism can allow us to 
go beyond Coxeter lattices and model other lattice uni¬ 
verses, something which would not be possible with the 
LW formalism because spherical cells, which that for¬ 
malism uses, can only well-approximate lattice cells that 
are regular polytopes [lsj . Once a particular arrange¬ 
ment of masses has been chosen, then by symmetry, the 
masses will maintain that arrangement on all subsequent 
Cauchy surfaces. The masses are therefore co-moving 
with respect to the Cauchy surface. 

Let us consider masses located in each tetrahedron at 


where Sli denotes the difference 5li := /,; + i — Zj. In terms 
of the strut-length nn, s 2 can be expressed as 


*? = 



(16) 


where we have made use of m to substitute for 5t 2 . 
Then varying Si with respect to nij yields 


dsj _ rrij 
drrij Si ’ 


(17) 


The area of any trapezoidal hinge between and 
is 


^trap 



+ h) 




(18) 


and varying this with respect to rrij yields 


dA T p t „ 

—7 - — — 77 TO i(‘i+1 + U) 

drrii 2 


-Xk+1-h) 2 - 


5ij. 

(19) 

Because all 4-blocks meeting at a trapezoidal hinge are 
identical, the hinge’s deficit angle simplifies to 


A trap = 2 tt - nOi 


( 20 ) 


where n is the number of faces meeting at the hinge and 
6i is the dihedral angle between any two adjacent faces. 
Since each trapezoidal hinge corresponds to the world- 
sheet of a tetrahedral edge and each face on this hinge 
to the world-tube of a triangle at this edge, n is equal to 










the number of triangles meeting at an edge; this number 
is given by the last column of Table HI 

To determine 9i , let us consider the representative 
hinge ABA'B'. Faces ABCA'B’C' and ABDA'B’D’ 
meet at this hinge and will be separated by a dihedral 
angle of 9i\ thus, we can determine Qi from the scalar 
product of the two faces’ unit normals. Let fi\ denote 
the unit normal pointing into ABCA'B'C' and fi 2 the 
unit normal out of ABDA'B'D'; then in co-ordinate sys¬ 
tem m, they have co-ordinates 


and 


(o.Q, 1 . 1 tv t) 
(i-Af 2 ) 4 




and their scalar product leads to the relation 


cos Oi = 


1 + H 2 

3-U 2 ' 


( 21 ) 


( 22 ) 


(23) 


We now have all the relevant geometric quantities nec¬ 
essary to solve the Regge equation (ll4l) . Thus by substi¬ 
tuting era d. m, cm, and (ED into ED. we are led to 
the Hamiltonian constraint 


li = 8 t tM ^ 

8 Ni 


H 2 - l 

8 1 

v 2 i 2 - 1 


1 


27 t — n6i 


(24) 


8 

Relation (l28l) can then be inverted to parametrise l in 
terms of 9, thus yielding 


l 2 = 8 


1 — 2 tan 2 



(29) 


Finally, the constraint equation (l25l) becomes 


l = 8 t tM 


N p tan | 

Nl [ 8?; 2 tan 2 § - ■§ (8v 2 


1 



2 - 7 T — n9 


(30) 


Equations (l29l) and (13011 provide a parametric description 
of the universe’s evolution in terms of 9. 

It is often easiest to study a model’s evolution by con¬ 
sidering the evolution of its Cauchy surface’s volume, and 
we shall do this below. In the continuum time limit, the 
CW Cauchy surface volume is given by 

TVq 

u N 3 {t) = -^=l(t)\ (31) 


where JV 3 is the number of tetrahedra in the surface and 
is given by the first column of Table U the volume’s rate 
of expansion is given by 


UN a (t) 


N 3 

2V2 


i(t) 2 i(t). 


(32) 


Both l and the lengths of the struts place constraints 
on the range of 9. In the continuum time limit, the 
strut-length is given by relation (1271) . and for the strut 
to remain time-like, we require that 9 > On the 
other hand, for l 2 to be positive, we require that 9 < 
2 arctan-^. Thus, 9 is constrained to lie in the range 


where Ni is the total number of tetrahedral edges on a 
Cauchy surface. By using (E3l) . we can express U as a 
function of the dihedral angle 9i and thereby re-express 

CUD as 


k = 8nM —^---2---. 

Nl [ 8 u 2 tan 2 | - i (8v 2 - 1 )] 5 “ n6 ' 

(25) 

We shall now determine the continuum time limit of 
the constraint equation. When taking this limit, such 
that as Sti —> 0 , 


9i —> 9 + O(dt ), 
li —> l, 

li+i —^ l T Idt -f- 0 (dt 2 ), 
the various geometric quantities become 
li — y 1 T O(dt), 

m 2 —> ^^' 2 — dt 2 + O^dt 3 ) , 

1 + I ; 2 

cos0j cos 9 ~-f-r—b O(dt). 

3-I / 2 


(26) 

(27) 

(28) 


TT 1 

— < 9 < 2 arctan ——. 
3 V2 


(33) 


For the square root in (l30l) to be real, we also require 
that 


9 > 2 arctan 



(34) 


2v 2 


We can com- 


thus l diverges if 9 = 2 arctan 
pare this with the range of 9, as given by (1331) . to 
see under what conditions will divergence occur. If 


2 arctan 


happens if 


2v 2 


> §, then l will definitely diverge; this 


Equality would correspond to placing the masses at the 
tetrahedral vertices. However the lower bound of 9 > ^ 
came from requiring that the struts be time-like, not from 
any direct constraints on l itself. Thus, even if the model 
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satisfies 2 arctan 



< |, we may still see the be¬ 


ginnings of a divergence in l that is abruptly cut off by 
the 9 = j bound. Therefore, for l to be unconditionally 
convergent, we require the more stringent constraint that 


0 <-i(8u 2 -l), 


which comes from requiring the square-root in (l30l) to be 
real and non-zero, even when 9 = 0; thus v 2 must satisfy 


When v 2 = 1/8, the masses are located at the mid-points 
of the tetrahedral edges. Therefore, the universe will only 
be unconditionally convergent if the masses are placed 
within a spherical region in the centre of the tetrahedron 
with a boundary that just touches the mid-points of the 
tetrahedral edges. Such a region would include the cen¬ 
tres of the triangles and of the tetrahedra itself. 

Interestingly, Collins and Williams 0 found a similar 
result in their study of closed dust-filled universes. They 
were considering different ways to measure the ‘time’ of 
the universe by using the proper time r of test particles 
located at different positions in the tetrahedron. We note 
that r is actually identical to the continuum time limit of 
Si, that is, the continuum time limit of the square-root 
of m . They found that dUN 3 /dr, where Un 3 is the 
volume of the universe as given by m, would diverge if 
the test particle were outside the same spherical region 
of convergence as the one we obtained. We suspect this 
region of convergence may be a generic feature of any 
model based on CW skeletons. 

Although masses situated at or near vertices will cause 
the resulting model to diverge, there may be a way 
around this problem. Each closed Coxeter lattice admits 
a dual Coxeter lattice centred on the original lattice’s 
vertices, as noted in Appendix [A] although the dual lat¬ 
tice may not necessarily be a tetrahedral lattice. Each 
model thus admits a dual model using Cauchy surfaces 
based on the dual lattice, and we can always extend the 
CW formalism to perform Regge calculus with these non- 
tetrahedral models. Where masses would have been at 
the vertices in the original model, when translated to 
the dual model, they would now be at the centres of the 
dual cells. As a result, the particles would now be both 
co-moving with respect to the Cauchy surface and follow¬ 
ing geodesics globally across the entire Regge space-time; 
this follows because even with non-tetrahedral cells, each 
cell would simply expand or contract uniformly about 
its centre as the Cauchy surface evolves, so a particle at 
the cell centre would simply propagate along the tem¬ 
poral direction only. We additionally conjecture that in 
non-tetrahedral models, there would also be an analogous 
spherical region of convergence centred at the cell centres 
regarding the placement of massive particles. However, it 
remains to be seen whether the spherical regions of both 
the original and dual models would completely cover the 


entire Cauchy surface such that for any configuration of 
the particles, there is always a model for which the con¬ 
figuration would be well-behaved. As noted in Appendix 
m the 5-tetrahedra lattice is dual to itself, so in this 
case, the region of convergence in the dual model is ex¬ 
actly identical to the region of convergence in the origi¬ 
nal model, and if masses are positioned at the vertices of 
the original model, then when they are translated to the 
dual model, the resulting model will behave in exactly 
the same way as the original model with the masses at 
its cell centres. 

In the foregoing discussion, we have only considered 
universes consisting of a single set of regularly distributed 
masses. We can also consider the more general case where 
we have not one but several different sets of masses, each 
set being regularly distributed and having equal mag¬ 
nitude, though the magnitude can differ between sets. 
If we work in the original lattice, then each set would 
be arranged in a way that preserves the lattice symme¬ 
tries, since, as mentioned above, the sets would lie at 
such sites as the lattice vertices or the mid-points of the 
edges. Therefore, even if the magnitudes of the masses 
differ between sets, because the lattice symmetries re¬ 
main preserved, the edges in particular would all remain 
identical to each other; hence, the lattice would still have 
only one length-scale l. A similar argument applies when 
the masses are translated to the dual model. Therefore 
for such models with multiple sets of masses, it can be 
shown that the continuum time equation for l becomes 

x—^ *, A/ tan it 1 

l =8 tt> Mi —A-3- -- 

i Nl [8^tan 2 |-i(8^ 2 -l)] 527r - n0 

(35) 

where i labels the set, the summation is over all sets, 
is the number of particles in set i, Mi is the magnitude 
of the masses in set i, and Vi denotes the parameter v 
for set i. It is clear that this equation has the same con¬ 
vergence conditions as those for the single-set universe. 
We note though that there are certain set combinations 
that would not lead to well-behaved models, even when 
translated to dual models. One example would be a 5- 
tetrahedra model where one set lies at the tetrahedral 
centres and another at the vertices; one of the two sets 
would always lie outside the region of convergence, re¬ 
gardless of whether the original model or the dual was 
being considered. 

Returning to the single-set universes, to illustrate the 
behaviour of our various Regge models, we have plotted 
dU/dt against U in Figure [3l where U is the universe’s 
volume as given by m for the Regge models. For com¬ 
parison, we have also plotted the corresponding graphs 
for the dust-filled closed FLRW model and for the equiva¬ 
lent LW models in 0 . Since the FLRW Cauchy surfaces 
are 3-spheres, the FLRW universe’s volume is given by 

^FLRw(i) = 27r 2 a(f) 3 , (36) 

and its rate of expansion by 

^FLRw(i) = 6 t r 2 a{t) 2 a(t); 


(37) 
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U(t) 


dust-filled FLRW . 600 tetrahedra Regge model - 600 tetrahedra LW model 

5 tetrahedra Regge model - 5 tetrahedra LW model 

16 tetrahedra Regge model - 16 tetrahedra LW model 

(a) 



U(t) 


dust-filled FLRW . 600 tetrahedra Regge model - 600 tetrahedra LW model 

5 tetrahedra Regge model - 5 tetrahedra LW model 

16 tetrahedra Regge model - 16 tetrahedra LW model 


(b) 

FIG. 3. The expansion rate of the universe’s volume dU/dt versus the volume U itself for the dust-filled FLRW universe, the 
three different Regge models, and the three equivalent LW models. In the Regge models, the masses have been positioned (a) 
at the centres of the tetrahedra and (b) at the mid-points of the tetrahedral edges. The Regge universe volume is given by the 
sum of the volumes of the Cauchy surface’s constituent tetrahedra, while the LW universe volume is given by the volume of a 
3-sphere with a radius equal to the lattice universe scale factor a(r). The universe’s total mass is the same across all universes. 


for a dust-filled FLRW universe, the scale factor a(t) is 
given parametrically by 


where 


a ° 1 1 \ 
a = y(l - COST?), 

t = y(?y- sin 77 ), 


a 0 


871-po 

3 


(38) 

(39) 


(40) 


and where po is the energy density when a = 1. We have 
taken the LW universe’s ‘volume’ to be the volume of 


a 3-sphere as well, though with the radius set equal to 
the LW lattice universe scale factor a(r) instead^] that 
is, the universe’s volume and expansion rate are given, 


2 We do not imply anything physical by this LW 3-sphere; it has 
merely been defined to facilitate comparison of length-scales be¬ 
tween LW and Regge models. We could instead have chosen to 
graph a(r) against a(r) for the LW universes and l/Z against 
l/Z, where Z is some suitably chosen constant that converts the 
Regge edge-lengths into an embedding 3-sphere radius; possible 
embeddings were explored at great length in [23 . The compari¬ 
son of the models would effectively be equivalent. 
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respectively, by 


Ulw(t) = 2t r 2 a(r) 3 , (41) 

Uvn(t) = 67 r 2 a(r ) 2 d(r), (42) 

and the LW scale factor a(r) is given parametrically by 


where 



r v„ r E b n , 

a = I —— (1 — cosrj), 



(V ~ sin ??), 


2 m 

rE ” = 

Eb = cos 2 V , 7 v 3 


(43) 

(44) 


and where the angle ^jv 3 is given by 

1 2 , 0jv 3 — sin 2'0/^ 3 

A 3 = 27T ’ 

with A 3 being the number of tetrahedra in the model, as 
listed in the first column of Table HJ Two classes of Regge 
models are shown in Figured on the one extreme, Figure 
I3al shows the behaviour when particles are positioned at 
the centres of the tetrahedra, while on the other, Figure 
I3bl shows the behaviour when they are at the mid-points 
of the tetrahedral edges, right on the boundary of the 
region of unconditional convergence. 

There are several features to note from these graphs. 
The evolution of all lattice models are well-behaved and 
closed. Moreover, the behaviour of the approximations 
converge to that of the FLRW universe as the number of 
particles is increased, which is to be expected because the 
matter content then becomes more dust-like. However, 
the Regge models and the LW models converge from op¬ 
posite directions. It is difficult to say a priori whether 
the LW approximation or the Regge approximation is 
more representative of the actual lattice universe: on the 
one hand, the average matter density is reduced as the 
number of particles is reduced, and this would weaken 
the gravitational interaction between particles, but on 
the other, the mass of each particle gets increased, and 
this would strengthen the interaction. However, compar¬ 
ison with exact initial value data on the time-sjunmetric 
hypersurface seems to favour the LW models [9(: using 
various measures for the lattice scale factor ai a tt on such 
a hypersurface, Clifton et al. have shown that the ratio 
O’latt /clflrw would always approach unity from above as 
the number of lattice masses increases. 

We note that there is a degree of ambiguity in choos¬ 
ing a common measure with which to compare the two 
approximations. It was demonstrated in [24J that CW 
Cauchy surfaces can be embedded into 3-spheres in E 4 ; 
the radius of the embedding 3-sphere could serve as a 
measure of the CW Cauchy surface’s scale factor, but it 
was shown that there were multiple ways to define the 


embedding and hence the radius, although each possible 
radius would simply be a constant re-scaling of the tetra¬ 
hedral edge-length l(t). Each radius implies a different 
3-sphere volume, which results in volume graphs of dif¬ 
ferent magnitudes, although regardless of which scaling 
is chosen, the graphs would still remain closed and well- 
behaved. Thus, it is within the realm of possibility that 
the models are actually equivalent and that the apparent 
difference is due to an inappropriate choice of scaling. 

However, we suspect that the principal factor behind 
the two approximations’ discrepancy is the finite resolu¬ 
tion of the lattice cell’s geometry in the Regge models. 
We have constructed Regge models where the entire ge¬ 
ometry of each cell has been reduced to one quantity, that 
of the tetrahedral edge-length; the cell’s interior geome¬ 
try consists uniformly of Minkowski space-time through¬ 
out, with no hinges or extra edges to provide addi¬ 
tional geometric information, and such a model is clearly 
very coarse-grained. In contrast, LW models approx¬ 
imate each cell’s interior geometry with Schwarzschild 
space-time, and, although not perfect, this approxima¬ 
tion should still be more accurate than completely uni¬ 
form Minkowski space-time. One could perhaps fine- 
grain the Regge models by further triangulating each 
lattice cell into smaller tetrahedra; this would certainly 
introduce more edges and hence more geometric infor¬ 
mation into each cell, although the smaller tetrahedra 
would no longer be identical nor equilateral. One algo¬ 
rithm for subdividing these tetrahedra has been devel¬ 
oped by Brewin |23j. Alternatively, one might be able to 
construct a finer-grained 5-cell universe, for instance, by 
distributing five massive particles in a regular manner in 
a 600-tetrahedra Cauchy surface; it remains to be seen 
whether such a regular distribution of the masses can 
be obtained. Though assuming it is possible, then each 
cell would now be associated with several edge-lengths 
rather than just one, and this would also provide greater 
geometric information about each lattice cell. Neverthe¬ 
less, we see that even these coarse-grained models can 
still reveal much qualitative information about the lat¬ 
tice universe’s behaviour, most notably its stability and 
its closed evolution; thus we shall continue to use them 
to study the lattice universe. 

Our first method for increasing the number of edge- 
lengths per lattice universe effectively ‘spaces-out’ the 
masses on the CW Cauchy surface so that each lattice 
universe cell would consist of multiple Cauchy surface 
cells. It may actually be easier to do this in Regge mod¬ 
els of flat and open lattice universes. The CW formalism 
should readily be extendible to Cauchy surfaces based on 
the flat and open Coxeter lattices described in Appendix 
[AJ Because such Cauchy surfaces would extend indefi¬ 
nitely outwards in all directions, it would be far easier to 
‘space-out’ the masses: for instance, in Cauchy surfaces 
based on the flat cubic Coxeter lattice, one could easily 
place a particle at the centre of every n th cube for some 
suitably large n, and each lattice universe cell would then 
consist of n 3 Cauchy surface cubes. In extending the CW 
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formalism to these other Coxeter lattices, it would also 
be interesting to see what analogous regions of conver¬ 
gence and dual models these new models would admit. 
For the moment though, we leave such investigations to 
future work. 

We shall henceforth work with universes where the 
masses are located at the centres of the tetrahedra. This 
would place the masses right in the middle of the region 
of convergence. Furthermore, as mentioned above, this 
is the only position in the tetrahedra where co-moving 
masses would also be following geodesics of the Regge 
space-time, which is closest to what happens in an FLRW 
universe. 


IV. PERTURBATION OF A SINGLE MASS 

We shall now construct a universe where the magni¬ 
tude of a single central mass gets perturbed from M to 
M'. This would induce a commensurate perturbation 
in the surrounding geometry, specifically in the lengths 
of the edges surrounding the perturbed mass. In gen¬ 
eral, there will be several sets of edges, each with its 
own independent length, and we must determine what 
those sets are. We begin with the edges in the tetrahe¬ 
dron enclosing the perturbed mass. The Cauchy surface 
will remain symmetric about this single perturbed mass; 
therefore the enclosing tetrahedron will still remain equi¬ 
lateral but with perturbed edge-lengths We shall 

refer to this tetrahedron as a Type I tetrahedron. Next, 
we consider the four tetrahedra that share a face with 
the Type I tetrahedron. We shall refer to these as Type 
II tetrahedra. These tetrahedra will no longer be equilat¬ 
eral: they will each have an equilateral base in the face 
shared with the Type I tetrahedron, but they will also 
have three other edges identical to each other but differ¬ 
ent from the l edges. In the 5-tetrahedra model, these 
two sets of edges are the only ones present, and we denote 
the length of the second set by /.■ 0) . In the 16 and 600- 
tetrahredra models, we must keep going; we must next 
consider the new sets of tetrahedra neighbouring the ones 
already considered; each additional set may introduce a 
new set of edges; and each additional set of edges will gen¬ 
erally have a length different from the lengths of all other 
edges so far considered. Thus, there will be far more 
than just two sets of edges in these models. For example, 
the 16-tetrahedra model has four distinct edge-lengths. 
One might consider simplifying the model by constrain¬ 
ing all edges to have a common length apart from 
the edges of the equilateral tetrahedron enclosing the 
perturbed mass; that is, one might attempt to localise all 
perturbation in geometry to just this equilateral tetrahe¬ 
dron. We attempted this but found that the resulting 
model was self-inconsistent. Therefore, we cannot study 
all three models in a general manner but must construct 
and examine each one individually. We have chosen to 
focus only on the 5-tetrahedra model as it is the most 


tractable, involving only two distinct edge-lengths. 

We shall once again work with a co-ordinate system for 
the tetrahedra’s 4-blocks, though each tetrahedral type 
requires its own system. Since the Type I tetrahedron is 
equilateral, we can use the same co-ordinates for its 4- 
block as those in (fTOl) and its E,; + i counterpart, although 
we must replace U and U+\ with and Z respec¬ 
tively; we shall continue to label the lower tetrahedron’s 
vertices by A, R, C, D and their upper counterparts by 
A', R', C', D'. 

For a Type II tetrahedron’s 4-block, we shall label the 
vertices of the lower tetrahedron by A, B , C, E, with 
E denoting the tetrahedron’s apex, and we shall denote 
their counterparts in the upper tetrahedron by A', B ', 
C', E'-, vertices A, B , C are shared with the Type I 
tetrahedron in E, and A', R', C' with the Type I tetra¬ 
hedron in Ej +1 . We assign to the lower tetrahedron the 
co-ordinates 


A 

R 

C 

E 



(i) 


2^3’ 

z. (1) 


2V3 


, 0 , 



°-§'°4 


(0,0,hi,0), 


where we have introduced the symbol 



(45) 


to denote the length of the tetrahedron’s central axis. 

By symmetry, the base A'B'C' of the upper tetrahe¬ 
dron should simply be a uniform expansion or contraction 
of the lower tetrahedron’s base; there may also be a shift 
dzi in the spatial direction orthogonal to the base, but 
that is all. This shift would be determined by the lengths 
of the struts. Thus, the upper base vertices A', R', C' 
should have the co-ordinates 


A' = 

B' = 

C' = 


' ,(i) /(i) \ 

Jj ±i _A±i Sz- i6t {1) 

2 ’ 2V3’ 1 y 

A(i) /(i) \ 

l J ±1 _A±1 Sz- i5t (1) 
y 2 ’ 2^’ ’’ 1 )’ 

(i) 


o, lj ±=,5zi, iSty 


where St ^ is the temporal shift of the upper base with re¬ 
spect to the lower one. Both Sz, and St^ can be deduced 
from the fact that (i) struts AA', RR', CC' are shared 
with the Type I tetrahedron and must therefore have the 
same length in both co-ordinate systems; and (ii) diago¬ 
nals such as AR', AC", BC' are shared with the Type I 
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tetrahedron and must therefore have the same length in 
both co-ordinate systems. From these two constraints, 
it can be deduced that Stf^ = Sti = U +1 — U and that 
5zi = but in the limit where M' -+ M, we 

must recover the unperturbed model, so we take Szi to 


be 6zi = ~lk Sl 


(i) 


The only constraint on apex E' is that it lie a dis¬ 
tance if' 1 from each of A', B', and C'. This is equivalent 
to the two constraints that E' lie at distance hi+i from 
the base’s centre and that the central axis connecting 
E' to the base’s centre lie orthogonally to the base. As 
base A'B'C' defines a 2-dimensional plane in a (3+1)- 
dimensional Minkowski space-time, the subspace orthog¬ 
onal to it would be a (l+l)-dimensional plane, and the 
tetrahedron’s central axis can be oriented along any di¬ 
rection in this plane. Combined with the first constraint, 
the second constraint implies that E' will lie on a hy¬ 
perbola in this (l+l)-dimensional plane; exactly where 
on the hyperbola it lies depends on the length of strut 
EE', which need not be identical to the length of AA!. 
Thus the axis of the upper tetrahedron may in general be 
Lorentz-boosted relative to that of the lower tetrahedron. 

Therefore, the upper vertices’ co-ordinates are given 
most generally by 



7(1) 


W X/(!) 

!*±I rgf 

2 ’ 2 ^ 3 ’ 2VG 


7« ,(1) e,(l) 

L i +1 b+1 0l j 

2 ’ 2 v / 3’ 2 1/6 


, iSti 


l 


( 1 ) 


lit! ,st 

V3 7 2^6’ 


61 


(i) 


E' = 0, 0, h i+ 1 cosh'i/’j — 


si l 1 ) \ 

^=, iSti + thi+i sinh xpi I, 


(46) 

where xpi is the relative boost between the upper and 
lower tetrahedra’s axes. 

However, as with the unperturbed model, we shall re¬ 
quire that all 4-blocks in the perturbed model have no 
twist or shear as well. This implies that the world-sheet 
generated by each tetrahedral edge between E, and E,; + i 
must be flat; in other words, the four vectors parallel 
to the four sides of this world-sheet must be co-planar. 
When this requirement is imposed on any world-sheet in¬ 
volving E', such as AEA'E', we obtain the further con¬ 
straint that 


0 = A^ — B^ cosh xpi + sinh^j, (47) 


This constraint effectively determines ipi, and solving it 
leads to a quadratic equation with two solutions. How¬ 
ever, if we take the limit where if' 1 -+ if ^ and 6zt —> 
— 51 f\ we expect to recover the 4-block of an equi¬ 
lateral tetrahedron; in this limit, xpi should vanish, and 
this only happens for one of the solutions. Therefore, we 
require that 


ipi = In 



Op 


(49) 


We note that because of the non-zero boost parameter xpi, 
the time-like hinges of this skeleton will be quadrilateral 
but not necessarily trapezoidal. Indeed, the lengths of 
the struts would not, in general, be identical: the length 
of EE’ would be different. Therefore, we shall denote the 
areas and deficit angles of quadrilateral hinges by A quad 
and <5j quad instead. We also note that in the continuum 
time limit where St j —> 0 , ipi must become xp dt to leading 
order in dt; this is because the relative boost between the 
lower and upper tetrahedra must become infinitesimally 
small as the separation between Cauchy surfaces tends 
to zero, and therefore, the zeroth order term of xpi must 
be zero. 

As mentioned earlier, Brewin has likened a choice of 
strut-lengths in the CW formalism to a choice of lapse 
function in the ADM formalism and the choice of hav¬ 
ing no twist or shear in the 4-blocks to a choice of shift 
function. In the ADM formalism, both the lapse and 
shift functions can be freely chosen independently of each 
other throughout the Cauchy surface. Yet in this CW 
skeleton, our choice of xpi implies that if we freely spec¬ 
ify the common length of struts AA', BB', CC ', then 
the length of EE 1 would be completely determined; we 
do not have any freedom to specify its length separately. 
We therefore see that the choice of shift function here has 
constrained the freedom to choose the lapse function to 
just the freedom to choose a single strut-length. Thus the 
CW formalism is not completely analogous to the ADM 
formalism, at least for global models. 

The global Regge action for the perturbed lattice 
model can now be expressed as 

8^5 = ^ A quad <5 quad + ^ Afd^ - 87 r M, Sij , 

^ (quadrilaterals ^ f triangular S i £ {particles} 

hinges / hinges / j £ {4-blocks} 

(50) 

where Sij is the path-length for the mass of magnitude 
Mi in the 4-block labelled j; when varied with respect to 
the struts to*,, it yields the Regge equation 


where 


A _ 7. 7 W 

|_i 7 

Bif, = hi + 1 f ^ 

CV = h i+ i/f 3 


hi +1 — 


iT 

2y/G 


( 48 ) 


° = £ 

i 


r)4 quad 

r quad 

dm k 1 


i £ {particles} 
j £ {4-blocks} 


(51) 


where the first summation is still over all quadrilateral 
hinges. 
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A. Global and local Regge equations 


As mentioned in Section [II] the global Regge equations 
can, in certain cases, be related to the local equations 
through a chain rule, and we have chosen to use such a 
chain rule relationship to obtain the global Regge equa¬ 
tions for the perturbed lattice universe model. We shall 
therefore elaborate now on this relationship. 

As discussed above, local variation requires that each 
4-block of the skeleton be fully triangulated. We have 
already described in Section [TT] how the triangulation of 
each 4-block is to be carried out. To understand how 
the triangulation works more globally, we label the five 
vertices of the 5-tetrahedra CW Cauchy surface E* by A, 
B , C, D , E, with E corresponding to the common apex 
of all non-equilateral tetrahedra; each vertex in E i+1 is 
then labelled by the same letter as its counterpart in £» 
but with a prime. Diagonals are chosen so that the letter 
of a diagonal’s vertex in alphabetically precedes the 
letter of the diagonal’s vertex in E*+i; thus, it is actually 
the relative ordering of the labels rather than the labels 
themselves that matter. In the 5-tetrahedra model, this 
triangulation generates the 10 diagonals of AB', AC', 
AD', AE', BC', BD', BE', CD', CE', and DE'. This 
algorithm for choosing the diagonals ensures that each 
4-block gets properly triangulated into 4-simplices; and 
indeed, the triangulation is consistent with the 4-block 
triangulation described earlier. 

To relate the local and global Regge equations, it is 
necessary that the global and local Regge actions be iden¬ 
tical. In our model, the Regge action will consist of two 
components, one corresponding to 


1 

87T 


E Aidi 

i £ {hinges} 


and another to 


E Mi Sij. 

i £ {particles} 
j £ {blocks} 


We shall show that each of these components are identical 
in the two actions and hence that the two actions are 
themselves identical for the model we are considering. 

We begin with the first component. In the global ac¬ 
tion, there are two types of hinges: the quadrilateral 
time-like hinges generated by the world-sheets of the 
tetrahedral edges and the triangular space-like hinges 
corresponding to the triangles of the tetrahedra in a 
Cauchy surface. In the 4-block described by fi5l) and 
(l46l) . an example of the quadrilateral hinge would be 
ABA' B', and an example of the triangular hinge would 
be ABC. In the local action, there are three types of 
hinges. The diagonals split each quadrilateral hinge into 
a pair of time-like triangular hinges, and these pairs of 
tinre-like triangular hinges correspond to the first type; 
an example pair would be ABB' and A A' B'. The second 
type are same triangular space-like hinges as in the global 


action. The third type correspond to isosceles triangular 
hinges formed by two diagonals and one tetrahedral edge, 
an example being ABC. These hinges have no analogue 
in the global action; however, they do not depend on 
the struts either, and as we shall be varying the action 
with respect to the struts only, they will have no con¬ 
tribution to the Regge equations; hence we shall ignore 
these hinges. Since each quadrilateral hinge is required 
to be planar, whether triangulated or not, the area of 
each quadrilateral hinge will equal the sum of the areas 
of its constituent tinre-like triangular hinges. Naturally, 
the space-like triangular hinges will have the same areas 
in the global and the local actions. Thus the hinge ar¬ 
eas appearing in the two actions are identical, and if the 
corresponding deficit angles are identical, then the first 
component is identical in the two actions. 

After edge-lengths of the same type have been set 
equal, the deficit angles Sf' and <5® of any triangular time- 
like hinge do become identical to the deficit angle ^ . quad of 
the original quadrilateral hinge; that is, = <5 quad . 

Setting the edge-lengths equal makes the two triangular 
hinges be co-planar both with each other and with the 
original quadrilateral hinge; thus the 4-blocks meeting 
at the triangular hinges would be flat; the unit normals 
of the triangulated faces would be identical to the unit 
normals of the original faces; and the dihedral angles be¬ 
tween the triangulated faces would be identical to the di¬ 
hedral angles between the original faces. Since the num¬ 
ber of faces meeting at the triangulated hinge is the same 
as the number of faces at the original hinge, the deficit 
angles for the triangulated hinge and the original hinge 
are identical. 

The deficit angle of the space-like triangular hinges are 
also identical in the original and the triangulated skele¬ 
tons; this follows simply because the unit normals to the 
faces meeting at the hinge will not change because of 
the triangulation. Thus we can conclude that the first 
component of the local and global actions are identical. 

As for the second component, clearly the path through 
the 4-block should not depend on whether the 4-block has 
been triangulated or not. Hence the global and local form 
of this component should also be identical; and therefore 
the local and global actions are indeed completely iden¬ 
tical. 

Since the two actions are effectively identical, we can 
express the global Regge equation as a linear combination 
of local equations through a chain rule 



dS dm\ 
^ dm f dmi 

3 1 



dS ddj 

ddj dmi ’ 


(52) 


where dS/dm\ is a local variation with respect to a single 
strut mf, with superscript l indicating that these struts 
are to be considered as local rather than global struts, 
and where dS/ddj is a local variation with respect to a 
single diagonal dj ; the first summation is over all local 
struts and the second is over all diagonals. 

We shall examine the contribution of each summation 
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in (Eg) to the global Regge equation. First, however, 
we note that as we are varying locally with respect to 
the struts, there are two distinct sets of struts in our 
skeleton relevant to local variation, one corresponding to 
the struts of the Type I tetrahedron, and the other to 
the EE' strut of the Type II tetrahedron. The first set 
of struts all have the same length, which we shall denote 
by m AA since the set includes strut AA'. The length is 
given by m as well, but with 4 replaced by l] 1 , that 
is, by 



(53) 


The length of strut EE' will be denoted by mf E and is 
given by 



hi+\ cosh-iA — hi -=dZ 

2VG 

1 2 

Sti + hi. |_i sinh^ 


(i) 


1 2 


(54) 


We now examine the contribution of the first summa¬ 
tion in (l52l) . We shall always choose the global strut- 
length rm to equal either mf E or m AA ; without loss of 
generality, let us assume then that rrii = m AA . Then it 
trivially follows that dm AA /drrii = 5ij, which is clearly 
0 (1) to leading order in dt when the continuum time limit 
is taken. If we substitute for Sti in (ED using ED, it can 
also be shown that the leading order of dm EE /drrii will 
be at least 0(1) as well. We shall be explicitly calculat¬ 
ing dS / dm\ below, and at the end of the calculation, we 
shall see that the leading order of dS/dmj is also 0(1). 
Therefore unless the second summation in (l52l) has neg¬ 
ative leading order, the first summation will definitely 
contribute to the leading order of the global Regge equa¬ 
tion. 

We now turn to the second summation in (1521) . Our 
skeleton has two diagonals, one triangulating quadrilat¬ 
eral hinges generated by l^ 1 edges, and the other tri¬ 
angulating quadrilateral hinges generated by 1 <A> edges. 
The first set of hinges are the ones that involve vertices 
E and E' in the Type II tetrahedra’s 4-blocks, and an 
example of a diagonal on such a hinge would be AE'; 
these diagonals have length d AE given by 



!('( i> ) 2 -(,; o> ) 2 +(<*); 

2h i+ i cosh^j- 1=51^ 

v 6 


(55) 


An example of a diagonal on the second set of hinges 
would be AB', and these diagonals have length d AB 
given by 

(d AB ') 2 = z£U (1) + (m AA ') 2 . (56) 


It can be shown that 

dd AB ' _ dd AB ' _ m AA ' 
drm dm AA ' d AB ' 

In the continuum time limit, where l ^ we have 

that m AA /d AB —> m AA /1^\ which is an O(dt) term. 
Thus, dd AB /drm raises the leading order of the second 
summation in (15211 by one. It can also be shown that 
in the continuum time limit, the leading order term of 
dd AE /drm is at least O(dt) as well. Therefore ddj/dnn 
will be at least O(dt) to leading order for all diagonals. 

So unless the leading order of dS/ddj is negative, the 
leading order of the second summation in (1521) will be 
at least O(dt). Naturally, verifying that the leading or¬ 
der of dS/ddj is not negative requires a direct calcula¬ 
tion of dS/ddj. However, such a calculation is beyond 
the scope of this paper. Instead, we shall assume the 
conclusion. First, if dS/ddj had negative leading order, 
then the corresponding Regge equation would diverge in 
the continuum time limit; thus our model would break 
down. We are assuming this is not the case. Secondly, we 
have found many similarities between the Regge model 
we are studying here and the parent Regge models of the 
A-FLRW universe studied in (24|. In that paper, we sim¬ 
ilarly considered the relationship between global and lo¬ 
cal Regge equations through an essentially identical chain 
rule to ED- We found the leading order of ddj/drrii to 
be O(dt) as well for all diagonals. We also found the lead¬ 
ing order of dm^/drrii to be 0(1); this followed trivially 
because all struts between pairs of consecutive Cauchy 
surfaces had equal length, so dm\/drrii would be unity 
for struts between the same pair of surfaces. Finally, we 
found the leading order of dS/dm[ to be 0(1) as well. 
Since the two models have identical leading orders for 
three of the partial derivatives appearing (l52l) and its 
analogue in [241] . we suspect they would have identical 
leading order for the final partial derivative, dS/ddj. In 
[2^ . the order of this term was O(dt), so we suspect it 
would be the same here. 

If this is true, then the second summation in (1521) would 
have a higher leading order than the first summation and 
would therefore not contribute. As a result, the solutions 
to 0 = dS/dmf would by themselves satisfy the global 
Regge equation 0 = dS/drrii, much like the situation 
with the A-FLRW Regge models in [24| . 


B. Particle trajectories 

The final term of the Regge action (1501) determines the 
effect of the masses on the behaviour of the universe. 
This term depends on s-jj , the length of the trajectory 
followed by the mass Mi through the 4-block labelled j; 
thus to fully specify our Regge model, we must specify 
what trajectory the masses will follow. 

Ideally, we should like our masses to follow geodesics 
throughout the entire universe and also be co-moving 
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with respect to the Cauchy surfaces, as we expect this 
to be the situation in the continuum universe. As men¬ 
tioned earlier on, this is the situation in the perfectly 
smooth FLRW universe, where test particles co-moving 
with respect to constant-t Cauchy surfaces are also fol¬ 
lowing geodesics of the space-time. In the lattice uni¬ 
verse, the point-masses should similarly be co-moving 
with respect to the universe’s Cauchy surfaces so as to 
preserve the lattice symmetries - we do not expect the 
gravitational interactions between a symmetric distribu¬ 
tion of masses to give rise to an asymmetric motion of 
the masses. Yet these particles should also be follow¬ 
ing geodesics as well. However we have seen that in 
the unperturbed Regge model, it is not possible for the 
particle to be simultaneously co-moving and following 
geodesics across the entire space-time unless the particles 
are positioned at the centres of the equilateral tetrahe- 
dra. Even in CW approximations of the perfectly homo¬ 
geneous and isotropic FLRW space-times, co-moving test 
particles will not follow geodesics globally either unless 
the particles are at the centres of the tetrahedra. Con¬ 
sider, for example, test particles co-moving with respect 
to the centres of the triangles. These centres themselves 
trace out piecewise linear trajectories as the underlying 
CW skeleton is piecewise linear. Therefore, test parti¬ 
cles co-moving with these points will follow the same 
piecewise linear trajectories, deflecting every time they 
cross from one Cauchy surface into the next. Only tra¬ 
jectories traced out by the centres of the tetrahedra will 
have no deflection, since the tetrahedra expand or con¬ 
tract uniformly about their centres as they evolve, and 
this thereby leaves the tetrahedral centres spatially fixed 
with time. We note though that each linear segment of 
a piecewise linear trajectory will still be a local geodesic 
within the 4-block it traverses because straight-line seg¬ 
ments are always geodesics according to the Minkowski 
metric. If we take the continuum space-time limit of 
a CW approximation to an FLRW space-time, we ex¬ 
pect to recover the continuum FLRW space-time itself in 
which any co-moving particle will indeed follow geodesics 
as well. The reason only the tetrahedral centres follow 
global geodesics is because Cauchy surfaces of the CW 
skeleton are not perfectly isotropic and homogeneous; 
thus, not all points on the Cauchy surface have been ‘cre¬ 
ated equal’. If co-moving particles do not follow geodesics 
in CW approximations of the perfectly homogeneous and 
isotropic FLRW space-times, there is even less reason for 
them to follow geodesics in approximations of the lat¬ 
tice universe, both perturbed and unperturbed. Thus 
we shall only require the point-masses of the perturbed 
Regge lattice universe to be co-moving with respect to 
the Cauchy surfaces. But in the continuum space-time 
limit, we do hope that these co-moving particles will in¬ 
deed follow geodesics as well. 

As mentioned previously, we have chosen to work with 
a lattice universe where the masses would be co-moving 
with the centres of the tetrahedra when the universe is 
unperturbed. For the particle in the equilateral Type I 


tetrahedron, the trajectory is straightforward: by sym¬ 
metry, the particle should remain, for its entire trajec¬ 
tory, at the tetrahedron’s centre. For particles in Type 
II tetrahedra, it is less clear where the particles should be 
positioned because the tetrahedra are no longer equilat¬ 
eral. The only clear symmetry here is in the equilateral 
base. We can therefore say that a particle should remain 
above the centre of the equilateral base for the entirety 
of its trajectory. The issue lies in fixing the particle’s po¬ 
sition above the base. We know that the vertices them¬ 
selves should be co-moving with respect to the Cauchy 
surface, so we shall use them as reference points to ex¬ 
press the trajectory of the co-moving particle. The par¬ 
ticle’s position on Cauchy surface E^ can be expressed 
as 

Pi = a( A + B + C) + /3E, (57) 

where A,B,C,E are the position vectors of vertices 
A, B , C , E, respectively, and where a and /3 are yet to 
be determined constants. For the particle to be inside 
the tetrahedron, a and /3 must be non-negative and sat¬ 
isfy the constraint 3 a + /? = 1. For the particle to be 
co-moving, we require that its position p i+1 on Cauchy 
surface E, + i be given by (1571) as well, but with vectors 
A, B, C, E replaced by A 7 , B 7 , C 7 , E', respectively. In the 
4-block between Ej and Ej + i, the particle then prop¬ 
agates in a straight line from p t to p i+1 , and such a 
trajectory is considered co-moving with respect to the 
Cauchy surfaces. 

There is one situation where there is clearly a unique 
choice for a and /3. Should become equal to if 3 ^ at 
any moment, then the corresponding tetrahedron will be 
equilateral; we would then require the particle to lie at 
the tetrahedron’s centre, which means a and (3 must be ^ 
at this moment. Yet based on our definition of co-moving 
trajectories, a and /3 must be constant over the particle’s 
entire trajectory. Therefore, a and /3 must be j over 
the particle’s entire trajectory. However, we shall take 
a and /3 to be 4 for all particles, regardless of whether 
their tetrahedra become equilateral or not. Such a choice 
would place the particles at the tetrahedra’s centroids, 
which would in some sense generalise our requirement 
that the particles be at the tetrahedra’s centres. 

C. Geometric quantities for the Regge equation 

We now turn to deriving the geometric quantities rel¬ 
evant for the Regge equation. From (l5ll) . it is clear that 
we need to derive three types of quantities: the varied 
areas of the time-like hinges, dA i /dmk', the correspond¬ 
ing deficit angles 6i , or equivalently, the dihedral angle 
&i between neighbouring faces; and the varied lengths of 
the particles’ trajectories across 4-blocks, dsij/drrik- 

We shall be taking the continuum time limit of these 
quantities so that we can express the Regge equation in 
its continuum time form. Thus, we shall be needing the 
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continuum time form of the lengths and the boost pa¬ 
rameters; these are given by 


lar hinges’ areas in terms of their edge-lengths. Hinge Al 
has area 


4 0) ->■ Z (0) (t), 

ll%^l i0) + l (0) dt + 0(dt 2 ) , 

4+ ) i-t/ (1) +Z' (1) dt + 0(dt 2 ) , 
ipi ip dt + 0(dt 2 } , 
d AE ' -a Z (0) + 0{dt ), 
d AB ' Z (1) + O(dt), 

—)■ m EE ' dt + 0(dt 2 ) , 
m AA -¥ m AA dt + 0(dt 2 ) , 


where 


h 


rh 


aa! 


rh 


EE' 


ij; 


Z(°)/(°) - §Z«Z« 
h 



-±=i™ - 
2 V 6 



hi « 


/lZ« 


MW 

i(i) • 
2^6 



(58) 

(59) 

2 2 

l) , (60) 
(61) 


Af 1 = - 


b (df s ') 2 (mfj 


(62) 


,\ 4' 


(<<■>) +(<ir) + (™A) 

and, when varied with respect to m AA , yields 


m 


AA' 


dA A1 

dm AA7 ~ 8/0 1 


(t ™) 2 + (<"') 2 - (mA)' 


Hinge B1 has area 


(63) 


= 1 


(©) 2 (<r)%(«o 2 

6 (<d B ’) 2 (-W)' 


(e) 4 +(<0‘ + H“')‘ 


(64) 


and, when varied with respect to mf A , yields 


and where h = 



(ZW) denotes the contin¬ 


uum time limit of hi. 

As mentioned previously, our skeleton has two types 
of time-like hinges corresponding to the world-sheets of 
l - n ^ and l- 1 ) edges. We shall refer to the quadrilateral 
hinge generated by as a Type I hinge and its trian¬ 
gular components as hinges Al and Bl, counterparts to 
triangles A and B, respectively, in Figure [I] We can use 


ti+l 



FIG. 4. A diagonal di divides a time-like quadrilateral hinge 
into a lower and upper triangular hinge, labelled A and B, 
respectively. The hinge has been represented with the generic 
tetrahedral edge-lengths of £; and A+i; these would equal l 
and respectively, or l and Z)i\ according to the type 
of hinge being considered. 

Heron of Alexandria’s formula to express these triangu- 


dA 


Bl 


771. 


AA' 


dm AA ' 8 A? 1 


(«,)A 0O 2 - H“') : 


(65) 

It can be shown that in the continuum time limit, the 
varied areas become 


dA 


Al 


dA 


Bl 


z« 


dm AA ' dm AA ' 


rh AA 


ktr )'- 1 


( 66 ) 

that is, dA A1 /dm AA and dA B1 /dm AA are identical to 
at least leading order. 

We shall refer to the quadrilateral hinge generated by 
z| 0) as a Type II hinge and its triangular components as 
hinges A2 and B2, also counterparts to triangles A and 
B, respectively, in Figure [2] Hinge A2 has area 


A A2 = 


(if) 2 (<iy ') 2 + (il 0 , ) 2 (mr)’ 
+ (<“') 2 

-'(if) 1 +(<-')* + ( mP ')‘ 

and, when varied with respect to m EE , yields 


(67) 


TO, 


EE' 


dA A2 
dmF E ' ~ 8Af2 


[}f ) 2 + {dt E ') 2 - (mf E j 


( 68 ) 
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Hinge B2 has area 


and, when varied with respect to m AA , yields 



. (70) 

It can be shown that in the continuum time limit, the 
varied areas become 


dAf 2 _ mf_{\ 


dA A2 

dmf E ' 


dA A2 

dm EE ' 



1 

3 






and 


dAf 2 
dm AA ' 


8A B2 

dm AA ' 



1 

3 





Since the Cauchy surface is no longer composed solely 
of equilateral tetrahedra, we must use ([6]) to calculate the 
hinges’ deficit angles. Each 4-block will contribute one 
dihedral angle to each of its hinges. The Type I tetra¬ 
hedron’s 4-block will contribute the same dihedral angle 
to all of its hinges: the tetrahedron is equilateral, so all 
of its associated Al hinges are geometrically identical, as 
are all of its B1 hinges; moreover, when all edges are con¬ 
strained to be identical, each pair of Al and B1 hinges 
become co-planar, as mentioned previously, and the unit 
normals to the two faces meeting at Al become identi¬ 
cal to the unit normals to the two faces meeting at Bl; 
thus all dihedral angles become identical, and this angle 
0-°^ is given by (l23l) but with replaced by l\ 1 ] . In the 
continuum time limit, we have that 6^ —> 9^°\ and 9^ 

is then given by G51) but with l replaced by that is, 
0 ( 0 ) 

is given by 


1 + 1 (iW) 

cos =-—- — + 0(dt). 

3 " i ('<■>) 


( 73 ) 


2 


1 

2 


+ O(dt) 


(71) 



The 4-block of a Type II tetrahedron has four distinct 
types of hinges: a pair of Al and Bl hinges for each 1 <A> 
edge and a pair of A2 and B2 hinges for each edge. 
We found that all hinges of the same type had the same 
dihedral angle, even though they are not entirely identi¬ 
cal because of the way the Regge skeleton is triangulated. 
Additionally, we found that Al hinges and Bl hinges had 
the same dihedral angle to leading order in the contin¬ 
uum time limit. Therefore in the continuum time limit, 
the 4-block of a Type II tetrahedron will contribute three 
distinct dihedral angles, which we shall denote 6A\ 0( 2 \ 
and 9^ for Al, A2, and B2, respectively. Each dihedral 
angle has been calculated by taking the scalar product of 
the unit normals to the two faces meeting at the corre¬ 
sponding hinge. 

One example of an Al hinge is ABB ', which is shared 
with the Type I tetrahedron’s 4-block. The two faces 
meeting at this hinge, ABB'C' and ABB'E ', are sepa¬ 
rated by a dihedral angle of d,- 1 "*. To leading order in the 
continuum time limit, this angle is given by 
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cos 9^ = 


(l O' 1 ’) 2 - 3 )“ 11 



h 




1 

2 


+ 0(dt). 


(74) 


The B1 counterpart hinge is A A' B ', and as mentioned 
previously, we found the two faces meeting on this hinge, 
AA'B'C' and AA'B'E’, to be separated by the same di¬ 
hedral angle as well at lowest order in dt. 


An example of an A2 hinge is AEE'. The two faces 
meeting at this hinge, ABEE' and ACEE', are sepa- 
rated by a dihedral angle of 6) '. This angle is given 
by 


COS 0 + = T 


(‘f’) 2 - ! (tf’) 2 ] (•nr’Y ~ hi 


§1 C 1 ) 

cosh .ipi hi 2 ^/g 


1 2 


(iff-iOi 11 )’ 


(raff) 2 - ft? 


/« 


^b+l COSh Ipi hi 2\f§ 


and in the continuum time limit, this expression becomes 


cosf^ 2 -* = —— 
2 


( Z (0)) 2 i (Z(D) 2 ] (m EE, y ~h 2 


h — 


jW ' 

2 V 6 


(/(0)) 2 _ 1 (z(i )) 2 (m EE ') 2 — h 2 


h — 


icn_ 

2 V 6 


0(dt). 


(75) 


Finally, the B2 counterpart to the above A2 hinge 
is AA'E'. The two faces meeting here, AA'B'E' and 

I 


AA'C'E ', are separated by a dihedral angle of 6^\ This 
angle is given by 


cos 


6 > (3) = 


\ 4+1 4 (1) - h >+1 ( ^75 cosh ipi + sinh ipi 




(«) 2 -4(lffi) 2 


-1 


J 2 4+1 4 1} + h i +1 ( ^75 cosh ipi + sinh ipi 


and in the continuum time limit, this expression becomes 


cos 9^ 


2 

1 h 

|_6 2 V 6 

2 

f m 

r- 

(Z(D) 2 _ 2 (/(°)) 2 " 


('■“) 

2 

+ 1 

(Z(i)) 2 _ 4 ( Z (o)) 2 ' 

)( 

'■->) 

2 

- 1 

+ 4 

12 z(1) + 2 V 6 h 

2 

O' 1 ’) 2 


+ 0(dt). 


(76) 


The last geometric quantity we require is the variation 
dsij/drrik of the particle’s path-length Sij through a 4- 
block with respect to each strut-length m*,. Between a 
pair of Cauchy surfaces 5+ and E? + i, there are only two 
distinct types of path-lengths: one corresponds to the 
path of the unperturbed masses, and we shall denote this 
path-length by s,;; the other corresponds to the path of 


the perturbed mass, and we shall denote this path-length 
by s'. It can be shown that varying s,; with respect to 
each of the relevant struts and then taking the continuum 
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time limit yields 


ds ds 
dm A dm B 

and 


ds 

dm c 


1 Tfl AA ( 1 ■ \ 

(77) 


ds 

dm E 


1 rh EE ' \hip + 1 
4 Sq hip + 1 


+ O(dt), 


(78) 


where m x denotes the length, in the continuum time 
limit, of the strut attached to lower vertex X , and where 
so denotes the quantity 


1 



This latter quantity is related to the continuum time limit 
of Si, where Si —>■ s as Sti —> 0, through the Taylor ex¬ 
pansion 


s ~ so dt + 0(dt 2 ) ; (80) 

the expansion has no zeroth order term because as 5ti —> 
0, the separation between Cauchy surfaces becomes in¬ 
finitesimal, and therefore the particle’s path-length from 
one surface to the next would become infinitesimal as 
well. It can also be shown that varying s[ with respect 
to each of the struts and then taking the continuum time 
limit yields 


ds' ds' ds' ds' i . 

dm A dm B dmP dm D 4 


(81) 


The derivation of these results has been explained at 
length in Appendix iBl 


D. Solving the Regge equations 


Having now determined all relevant geometric quanti¬ 
ties, we can now substitute them into (1511) to obtain the 
corresponding Regge equations. As we have two distinct 
types of struts, locally varying the Regge action will lead 
to two distinct equations, 


E 

i 

E 

i 


d %5 t = j2^ M > dSiJ 


dm EE ‘ 


dm EE '' 


8 \4 = 5 > m . dsv 


dm AA 




dm AA ' 


(82) 

(83) 


We can directly obtain the continuum time limit of these 
equations by substituting in the continuum-time form of 
the geometric quantities, and this is how we shall pro¬ 
ceed. 

We begin with the first equation. Between each pair 
of consecutive Cauchy surfaces, there is only one strut 


with length m EE , namely the strut EE' in the Type II 
tetrahedron’s 4-block. Thus the only relevant geometric 
quantities are those involving strut EE' and vertex E. We 
begin by working out the left-hand side of ( 1571 ) . starting 
with the quantity - —^* E / . The only edges meeting at ver- 

° m k 

tex E are AE- type edges, all of which have length 
Thus the only hinges meeting at strut EE' are hinges 
like AEE' , what we have called A2 hinges above, and for 
each of these hinges, -^ew is given by CD in the con¬ 
tinuum time limit. Next, we consider the corresponding 
deficit angle. In the 5-tetrahedra model, three faces meet 
at each hinge, and hence three dihedral angles contribute 
to the hinge’s deficit angle. The only dihedral angle at 
A2 hinges is 6^ 2 \ which is given by ( 1751 ) in the continuum 
time limit; thus the deficit angle is Si = 2 tt — 3 0( 2 \ We 
finally perform the summation on the left-hand side of 
( 1571 ) . Because Si is identical for all hinges meeting 

at EE', performing the summation is equivalent to mul¬ 
tiplying this term by the number of hinges at EE'. As 
there are four edges meeting at vertex E , there can only 
be four hinges. 

Next, we consider the right-hand side of (|82|) . There 
are four tetrahedra meeting at vertex E, and each carry 
an unperturbed mass M; thus Mi = M for all i. Since 
all four masses and all four tetrahedra are identical, the 
quantity drn E E > w ill be identical for all i as well; when 
index j = k, its continuum time form is given by CD); 

ds ■ ■ 

otherwise it is zero. Because Mi -—^ is identical for all 

dm k 

i , performing the summation on the right-hand side of 
(1521) is equivalent to multiplying M EE , by four. 

° m k 

Finally we substitute everything into ( 1571 ) to obtain the 
constraint equation 


Z (0 > (2ir- 3(9 (2) ) 


jhip + 1 

= AttM ———- 


So 


z«\ 2 

m) 


2^6 


- h 


hip + 1 


- 1 


(84) 


We next consider the second Regge equation (l83l) . 
Struts of length m AA are connected to the four vertices 
labelled A, B, C , D in the Type I tetrahedron. As we 
shall see, varying any of the four associated struts will 
lead to the same Regge equation. We begin with the left- 
hand side of (1551) as well. Each of the four vertices has 
four tetrahedral edges attached to it: three edges are at¬ 
tached to the three other vertices in the tetrahedron and 
have length l W; the fourth is attached to vertex E of a 
Type II tetrahedron and has length . Each of the four 
edges generates a pair of time-like triangular hinges, one 
of which is attached to the vertex’s strut. The length-^ 0 ) 
edge contributes hinges like AA'E' in a Type II tetrahe¬ 
dron’s 4-block, what we have called B2 hinges above, and 
dA 

for such hinges, a 'a> is given by (1721) . The other three 
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hinges correspond to either Al or B 1 hinges, and in ei¬ 
ther case, - cL \ i A , is given by (1M1) . Once again, the deficit 

angles at each hinge involve three dihedral angles. At a 
B2 hinge, the only relevant dihedral angle is which 
is given by GO). So a B2 hinge has a deficit angle of 
Si = 2ir — 3 6( 3 \ At Al and B 1 hinges, the relevant di¬ 
hedral angles are and 9^\ In general, each dihedral 
angle at a hinge comes from a 4-block meeting at the 
hinge: each 4-block has two faces meeting at the hinge, 
and the dihedral angle contributed by the 4-block would 


be the dihedral angle between these two faces. Addition¬ 
ally, each 4-block is generated by a tetrahedron attached 
to the edge that generates the hinge. So in the case of Al 
and B1 hinges, each length-^ 1 ) edge is always attached 
to the Type I tetrahedron and to two Type II tetrahe- 
dra; the 4-block of the Type I tetrahedron contributes 
a single 9^°' ) to the hinge’s deficit angle; the two Type 
II tetrahedra each contribute a $T) angle; so the deficit 
angle of both Al and B1 hinges is Si = 2n — 9^°' ) — 2 9^\ 
Combining the contributions from all four hinges, we can 
express the left-hand side of (1551) as 


E 


dA z 

dm AA 


7 Si = 


lA)m AA ‘ 



(2tt-(9 (0) -2 6 ) (1) ) 


+ 



1 

3 




2 




We now move on to the right-hand side of (1551) . Vary¬ 
ing the strut-length will affect the trajectory length sy 
of four neighbouring masses. One of these will be the 
perturbed mass of magnitude M' while the other three 
will be masses of magnitude M. Thus the right-hand side 


will be 


13 


= 87 r 


M' 


ds' 

dm A 


3 M 


ds 

dm A 


ds' 


with the quantity given by (1511) and -^x by (1771) . 
Finally substituting everything into (1551) . we obtain 



where the quantity l/sq appearing on the left-hand side 
would actually be real; this follows because if the path- 
length s of the unperturbed particles is time-like, then s 
would have to be imaginary, and by virtue of its relation¬ 
ship to so through Taylor expansion (15U1) . So would have 
to be imaginary as well. We note that both equations 
(1551) and (1551) came from the 0(1) term of the Regge 
equations 0 = g ^ E > and 0 = g ^f A , ; thus, the Regge 

equations are 0(1) to leading order, as claimed in Sec¬ 
tion IIV Al when we were relating the local and global 
Regge equations through chain rule (1571) . 


These Regge equations however involve both 1 *-°) and 
/7) i n a non-linear manner, which makes solving for them 
difficult. We shall therefore linearise these and all subse¬ 
quent equations by performing a perturbative expansion 
up to first order in SM := M'—M. Under this expansion, 
we must have that 

Z< 0 > ^l+Sl w , 

/ (0) rj i + si w , 

l w &l+SlW, 
i w asi + si w , 
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ip « Sip, 
ip w Sip, 

8® t* 0 + 50i for i = 0,1,2,3, 


where we have made use of (1771) and (1571) to express these 
as functions of 0. We note that only S9 2 depends upon 
the boost parameter ip; if we substitute in relation (1551) 
for Sip, then S9 2 becomes 


as the zeroth order terms must match the corresponding 
quantities for the unperturbed model. The zeroth order 
angle 9 is given by relation (1751) . 

It can then be shown that the zeroth order terms for 
both Regge equations yield 


47 tM 
2tt — 39 



( 86 ) 


which is equivalent to the unperturbed Regge equation 
m for the 5-tetrahedra model with the masses at the 
tetrahedral centres: for the 5-tetrahedra model, we would 
have N p = 5, N = 10, and n = 3; and since the masses 
are at the tetrahedral centres, we would have v 2 = 0, 
because recall that v is the ratio between \vi\, the dis¬ 
tance of a mass to its tetrahedron’s centre as given by 
(fl5l) . and k, the tetrahedral edge-length. By using (1771) 
to substitute for l , we can parametrise l entirely in terms 
of 9, yielding 


l = 


4\[2-kM 
2tt — 39 


tan 



(87) 


An expression for Sip can be deduced by taking the 
perturbative expansion of ip as given by (ED; we thus 
obtain 


Sip = 




( 88 ) 


cos 9 — Sl^ 1 ^ 

is9)Si {1) , 

(93) 

and we note that SI^ has now dropped out of this ex¬ 
pression; in fact, now none of the angle perturbations 
depends on Sl^°\ 

From the perturbative expansion of (1541) . the first order 
term yields 


S9 2 = - 


4*j2 


(2n — 3 9) (1 + cos 


irM 


1 — cos 


3 cos 9 — 1 
1 — cos 9 


(1 + 


Sl {0) (2tt - 3 9) - 31S9 2 


7 tM 1 j i 

It 


(94) 


where we have substituted for Sip using (1881) . and from 
the perturbative expansion of (1851) . the first order term 
yields 


4ttSM 



+ 3 

(2tt — 3 9) 

1 - A P 


(27t — 3 9) — 31 (S9 0 + 2S9 X + S9 3 ) 

U 2 (si {0) - si (1) ) + y i i si w . 

(95) 


The quantities S9i can be deduced from the pertur¬ 
bative expansions of relations (1731) to dZHD- The zeroth 
order terms of these relations are all identical to (1751) . as 
expected. The first order terms yield 


59 * = - 


59i = - 


Sft 1 ' 1 1 3 cos 9 — 1 


4V2 

si w 


1 — cos 9 
3 cos 9 — 1 


(1 + COS0) , 


(89) 


(1 + cos 9) 


4\f2 V 1 — cos 9 
1 (27t — 3 9) 

--7T- (1 T COS 9 ) 

8 ^2 nM y J 


(<5Z (0) - « (1) ) , (90) 


= ( 2 *-3*) ( (2cos9 _ !) 

4 -JIkM Vi-cos e) ■ 

,3cO8 ''“ 1 (l + cos9)(24i‘ 0 >- 


1 

4\/2 V 1 — cos 9 
4 7 tM 


S9 3 = - 


V3 (2tt-3 9) 

(27r — 3 9) /1 + cos 


(3cos0 — 1) Sip, 


4^2 ttM V 1_cos 

SiS 11 / 3 cos 9 — 1 


cos 9 


4V2 V 1 — cos ( 


(1 + cos 9), 


(91) 


(92) 


We note that none of these relations depends on SI^ 
either. 

In the unperturbed model, we used the dihedral angle 
9 to parametrise l through relation (1791) . We shall do 
something similar here and parametrise SI ^ with respect 
to one of the angle perturbations. It is easiest to do this 
with relation (1571) . which then yields 


« (1) 



1 — cos 9 S9 0 

3 cos 9 — 1 (1 + cos 9) 


(96) 


Since none of the angle perturbations depends on Sl^°\ 
a similar parametrisation is not possible for 5l(°\ How¬ 
ever, the first order terms (1771) and (1751) of the two Regge 
equations (1571) and (1551) do not depend on 51^ anyway, 
so such a parametrisation of SI^ is not necessary. 

Using relations (1571) . (1701) . (1751) . and (1771) to substitute 
for S9 0 , 59\, 59 2 , 59 3 , relation (1701) to substitute for 51^, 
and relations ED and (1771) to substitute for l and l, we 
can now solve (1771) and (1701) for 51^ and 51 ^ and express 
them exclusively in terms of the parameters 9 and 59 3 . 
We find that 
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Sl w = 


4\/2irM 
3(2tt-3 0) 


2sin0 | 1 + 2cos6 ) -f (271- — 3 0) 


1 + cos 0 


i -i 


SM 

1 — cos 0 

1 

2 

6 sin 0 

COS0 

M 

1 + COS0 


1 + cos 0 


+ (2tt-30) 


3 cos 0 — 1 


3 (50n 


and that 


Sl w = 


(27r — 3 9) 1 + cos 9 

4:V2ttM 


3sin0 + (27 t — 3 0) 


2 sin 9 


1 + cos t 

1 + 2 cos 9 
1 + cos 9 


3(2tt-30) 


2 8 i„f>| 1+2cos8 |+(2 T -3fl) 


1 + cos 9 


SM 

1 — cos 0 
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2 

6 sin 0 

COS0 

M 

1 + COS0 


1 + COS0 


+ (2tt- 3 0) 
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3(50n 


(27t — 3 0) 1 + cos 0 


3sin0 + (27 t — 3 0) 


2 sin 0 


1 + 2 cos < 
1 + cos 0 


(97) 


(27T-3 0) 


(98) 


(2tt- 3 0) 


Finally, we note that our two parameters 0 and S9q are 
not independent of each other; rather both are functions 
of the underlying time parameter t; therefore if one pa¬ 
rameter evolves, so must the other. We can relate the 
two parameters to t through the system of differential 
equations 


(5/ (1 ) = -f-(5Z (1) , (100) 

dt 

where the left-hand side denotes the quantities given by 
m and (19611 . while the right-hand side denotes the ex¬ 
plicit differentiation of (1571) and (1551) with respect to t. 
The first equation involves only 0; it can be solved on 
its own to yield 9(t). This can then be substituted into 
the second equation to give a differential equation for 59$. 
We shall solve these equations numerically. To determine 
a unique solution though, we must also specify a set of 
initial conditions: we shall require the perturbed model 
to obey the initial value equation at its moment of time 
symmetry. This equation will be explained in the next 
section, and it implies specific conditions on 0 and 59q 
which we shall derive. 

The range of the parameter t will be constrained by 
the requirement that all struts remain time-like. As this 
constraint depends on SM, the resulting range of t will 
also depend on SM. 

Before leaving this section, we wish to remark on an 
advantage that local variation has afforded over global 
variation. For each model, local variation has yielded a 
pair of Regge equations that, when expanded perturba- 
tively, gave three distinct equations, an identical equation 
from their zeroth order terms and two distinct equations 
from their first order terms. Had we directly varied the 
action globally instead, we would only have obtained one 
Regge equation corresponding to a linear combination 


of the two local Regge equations, and the perturbative 
expansion of this global equation would give just two in¬ 
dependent equations. Thus local variation has provided 
us with an extra independent equation, allowing us to 
specify one more of the five quantities, l, l, Sl^°\ Sl^\ 
Si (•*■), that we needed to solve for. This has allowed us 
to pick one unique solution out of many in the solution 
space of the global Regge equations. 

V. INITIAL VALUE EQUATION FOR 
PERTURBED MODELS 

We now have a set of equations, in (1551) and (11001) . that 
should determine the entire Regge skeleton from a set of 
initial conditions on 0 and 80q. However, the question 
remains as to what initial conditions would be appro¬ 
priate. To answer this, we shall consider the analogous 
(3+l)-formulation of general relativity wherein the entire 
space-time is similarly determined by evolving forwards 
in time from some Cauchy surface So with initial data; 
naturally, the evolution equation would be derived from 
the Einstein field equations. It has been shown |25j that, 
for such a formulation, the required initial data consists 
of the first and second fundamental forms, h and X'i the 
former corresponds to the projection of the metric g into 
Eo and effectively determines the 3-dimensional intrin¬ 
sic curvature of So; the latter effectively determines the 
extrinsic curvature of Eo within the overall space-time. 

However, for this initial data to be consistent with the 
Einstein field equations, there is a set of constraint equa¬ 
tions that it must satisfy. Let us express the Einstein 
field equations in the form 

G = 8ttT, 

where G is the Einstein tensor and T the stress-energy 
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tensor. Let n denote a field of normalised one-forms ev¬ 
erywhere orthogonal to E 0 . By making use of the Gauss 
equation 

^ ^ Rpuap RalS'yS h fj,tl v h? a h p Xpc t Xvp d - Xpp Xiscn 

( 101 ) 

which relates the 3-dimensional intrinsic curvature 
^Rpvap of So to its extrinsic curvature x and its 4- 
dimensional intrinsic curvature R a /3-yS, we can express 
the relation G (n, n) = 8nT (n, n) as 

Mr + (xT V ) 2 - XT X pa Kp h va = 167 TP, (102) 

where ^ R is the 3-dimensional Ricci scalar of Eo and P 
is the energy density of the matter source as measured 
by an observer co-moving with respect to So- Equation 
(11021) gives the first constraint equation; it is actually the 
Hamiltonian constraint of the ADM formalism, where it 
is customarily derived by extremising the ADM action 
with respect to the lapse function [26|. Let {«i}, for i = 
1,2,3, denote a set of normalised basis vectors tangent to 
Eo, and let | denote covariant differentiation with respect 
to the metric connection implied by h. By making use 
of the Gauss-Codazzi equation 

Rep n°h p M = X a p\a ~ X a a\pi ( 103 ) 

which relates the extrinsic curvature x of Eo to its 4- 
dimensional intrinsic curvature in the form of the Ricci 
tensor R ap , we can express the relation G(n, tij) = 
8nT (n, itj) as 

(x CT/ V ~ hap ) = 8?r Tpv n>1 U ”( 104 ) 

which is actually a set of three equations, one for each 
i. This gives the rest of the constraint equations. These 
are the momentum constraints of the ADM formalism, 
where they are customarily derived by extremising the 
ADM action with respect to the shift functions (26| . 

Quite often, the initial surface Eo is chosen to be the 
surface at a moment of time symmetry, that is, the mo¬ 
ment when the surface’s extrinsic curvature, as given 
by the second fundamental form Xi vanishes. In this 
case, the momentum constraints would vanish while the 
Hamiltonian constraint would simplify to 

(3) f? = 16?rp; (105) 

this is known as the initial value equation at the moment 
of time symmetry. 

To determine the appropriate initial conditions for our 
Regge model, we shall require our model to satisfy this 
equation as well at its moment of time symmetry. How¬ 
ever, we note that this equation and the Einstein held 
equations, from which it is derived, will only be satisfied 
in an average manner on a Regge Cauchy surface. Curva¬ 
ture in the surface is concentrated only at the hinges, yet 
matter can be distributed away from the hinges where 
the skeleton is flat, which is indeed the case for our 


Regge model; thus the two sides of the equation will 
not agree in a point-wise manner. This contradiction 
arises because the Einstein held equations actually apply 
to smooth manifolds rather than Regge skeletons; they 
come about by varying the Einstein-Hilbert action when 
the underlying manifold is smooth rather than discrete. 
Thus by using the Einstein equations in this manner, we 
are effectively varying the Einstein-Hilbert action on a 
smooth manifold hrst and then applying the resulting 
held equations on a discrete manifold afterwards. The 
standard approach in Regge calculus is to use a discrete 
manifold from the very beginning, with the held equa¬ 
tions obtained being different as a result. Clearly, the 
two approaches are not equivalent. 

We shall now apply (11051) to our Regge model and 
thereby deduce a set of initial conditions on 9 and 59q. 
The perturbed models attain a moment of time symme¬ 
try when all lengths cease expanding or contracting, that 
is, when = l T) = 0. We shall perform an averaging of 
(I105D by integrating it over such a time-symmetric CW 
Cauchy surface Eo- The left-hand side of (11051) then be¬ 
comes [2(1 

f {3) Rd 3 x = 2 

i f= {hinges} 

where the integration measure is unity because the Regge 
tetrahedra are hat, the summation is over all edges in Eo 
because these would be the hinges of a 3-dimensional 
skeleton, £, is the length of an edge, and Si is its corre¬ 
sponding 3-dimensional deficit angle; the right-hand side 
evaluates to 



Therefore, the averaged initial value equation for the per¬ 
turbed model can be expressed as 

kSi = 8tt (5 M + SM). (106) 

{edges} 

By requiring our model to satisfy this form of the ini¬ 
tial value equation, we shall deduce the necessary initial 
conditions on 9 and 59q. 

The only quantities in (11061) that have yet to be de¬ 
termined are the deficit angles. A Cauchy surface of the 
perturbed universe has only two distinct types of hinges, 
the edges of length and the edges of length l ^. Each 
edge is connected to three faces separating three tetrahe¬ 
dra, so each tetrahedron at the edge will contribute one 
dihedral angle to the edge’s deficit angle. A Cauchy sur¬ 
face of the perturbed universe also has only two distinct 
types of tetrahedra, the Type I and Type II tetrahedra. 
As the Type I tetrahedron is equilateral, it will contribute 
the same dihedral angle to each of its six edges, and we 
denote this angle by . In the Type II tetrahedron, all 
edges of are identical to each other, as are all edges 
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of length Z-'b Thus this tetrahedron will contribute the 
same dihedral angle <f)^ to each of its Z- 0 ^ edges and the 
same dihedral angle <^- 2) to each of its edges. To de- 

termine the dihedral angle cf)}' between any pair of faces, 
we shall again take the scalar product of the unit normals 
to the two faces, as this product will yield cos<^. 

Let us first consider the dihedral angles in the Type 
I tetrahedron. We can use co-ordinate system (flOl) for 
this tetrahedron, dropping the time co-ordinate so that 
we work in a purely 3-dimensional spatial co-ordinate sys¬ 
tem and replacing lengths U with l^\ as this tetrahedron 
has edges of length Z^; we can then use this co-ordinate 
system to assign co-ordinates to any of the normal vec¬ 
tors. We can calculate using the faces meeting at 
edge AB; these faces are ABC and ABD , and the scalar 
product of their unit normals yields 

cos# = i. (107) 

For the Type II tetrahedron, we can work with 
co-ordinate system (1451) . again dropping the time co¬ 
ordinate to obtain a purely 3-dimensional spatial system. 
We can use edge AE to calculate the dihedral angle <p\ 11 
at an edge of length Z^; the faces meeting at AE are 
ABE and ACE , and the scalar product of their unit 
normals yields 


, (l) 

cos q>\ ’ = 



(108) 


Similarly, we can use edge AB to calculate the dihedral 
angle ] of an l 1 - 1 ' 1 edge; the faces meeting at AB are 
ABC and ABE , and the scalar product of their unit 
normals yields 


To calculate the deficit angles at an edge, we simply 
subtract the three relevant dihedral angles from 2-7T. Since 
only Type II tetrahedra have Z^°l edges, an edge must 
be connected exclusively to Type II tetrahedra, and its 
deficit angle 5^ must therefore be 

= 27t — 3 (j)^ 

ft 27t — 3 arccos ^ (si^ — Sl^ 1 ^ . 

As there is only one Type I tetrahedron on the entire 
Cauchy surface, an ZW edge can only be connected to 
one Type I tetrahedron, and its other two tetrahedra 
must be Type II; its deficit angle must therefore be 

<5(°)=27r —</>(°)-2</>( 2 ) 

ft 27T - 3 arccos(« (0) - SI (1) ) . ^ 


Having now determined the deficit angles, we can now 
substitute all relevant geometric quantities into (11061) . As 
mentioned previously, a 5-tetrahedra Cauchy surface will 
have a total of 10 edges; six must come from the Type 
I tetrahedron and must therefore be of length Z^ 1 ); the 
remaining four must therefore be of length Thus 

(I106D can be expressed as 


8tt (5 M + SM) 
= 4 Z (0) 


„ o 1 V2 

2n — 3 arccos - H—— 

O L 


0 i 2V2 

27T — 3 arccos-— 

3 3 Z 


-6 Z (1) 

10Z + 2(2<5Z (o) + 3(5Z (1) ) 


2-7t — 3 arccos ■ 


1 

3, 
(114) 


/ ( 2 ) 

COS (j)\ 1 = 



(109) 


To obtain the continuum time limit of these expressions 
to leading order in dt , we simply need to drop all sub¬ 
scripts i. If we next take the perturbative expansion of 
these continuum time expressions, such that 

(j)^ « cj) + 

(j)^ ft! (j) + 

we find that 

<t>=\, (HO) 

and that 

S(j,A) = -$(/>( 2 ) = - ^ - Sl w ^j . (Ill) 


The zeroth order term corresponds to the initial value 
equation for the unperturbed model. At the moment of 
time symmetry in the unperturbed model, we have that 
1=0, and from (1281) . it follows that 8 is 


9 = arccos -. 

3 


(115) 


Substituting this into the unperturbed model’s Regge 
equation (15()1) . we deduce that 


47tM = Z 


27T — 3 arccos - 
o 


(116) 


which is identical to the initial value equation for the un¬ 
perturbed model. Thus the Regge equation of the unper¬ 
turbed model satisfies its initial value equation at the mo¬ 
ment of time symmetry. Therefore, for the zeroth order 
component of the perturbed models’ initial value equa¬ 
tion to be satisfied, we must also require that 8 satisfy 
condition (11151) . This is the initial condition on 8. 
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We can next deduce the condition on 59q by solving for 
it from the first order term of (llldl) . After using (11151) 
to substitute for 6 as well as (liiTl) and (j98)) to substitute 
for 51 and 5l^\ we find that 58o must satisfy 

56 o = 0. (117) 

We note that the behaviour of SI W near time symme¬ 
try cannot be determined from (l96l) . On the one hand, 
condition (I117D suggests it may approach zero, while on 
the other, condition (11151) suggests it may diverge. As 
we shall see below however, 51^ is indeed well-behaved 
and tends towards zero as time symmetry is approached. 

Finally, even after imposing conditions (11151) and (11171) 
at some moment t = T max , we must still ensure that 
51 = 51 W = 0 at that moment; otherwise, t = T max 
would not be a moment of time symmetry. We have just 
mentioned that 5l l[ ' 1 does tend toward zero as time sym¬ 
metry is approached. Given this information, it can also 
be shown that 51 will be zero as well. First, by com¬ 
paring m and (1551) . we note that 51^ can be expressed 
in terms of 51 ^ to give 



FIG. 5. The expansion rate of the universe’s volume dU/dt 
versus the volume U itself. Graphs corresponding to four 
different values of 5M/M are shown, with 5M/M = 0 cor¬ 
responding to the unperturbed model. Owing to 8M/M = 
0.001 and 5M/M = 0.01 being very small perturbations, the 
graph for <5 M/M = 0.01 has almost completely covered the 
graphs for 5M/M = 0 and 8M/M = 0.001. The left of the 
graphs have been truncated at the moment the struts turn 
null. In all four models, the mass M has been fixed to be 
M = 1/5. 


^( 0 ) _ ^(i) _ 8^7 t5M ' 1 - cos6> ~ 2 
3 1 + cos 6 

1 -l 

+ (2tt-3(9) . 

(118) 

If we differentiate this with respect to t, we obtain an 
expression of the form 

<5/(°) = fl'Ci) - 8 ^ n6M Fl (0) 0. 

3 

At the moment of time symmetry, we have said that 
51 W will vanish and that condition (11151) must be sat¬ 
isfied; then it can be shown that F 1 (9 = arccos|) will 
not vanish, and therefore 51^ will vanish if and only if 9 
vanishes. To see that 9 does indeed vanish, we next dif¬ 
ferentiate l, as given by (ED, with respect to t to obtain 
an expression of the form 

i = F 2 {8)8. 

The left-hand side is given by (l29l) and is zero at time 
symmetry. It can also be shown that F 2 {9 = arccos i) 
will not vanish; thus for the two sides of the expression 
to equal, it follows that 9 must vanish at time symmetry. 
Therefore, we deduce that 51^ vanishes at time symme¬ 
try provided 51^ vanishes as well. 



FIG. 6. The same graph as in Figure [o] for the 5M/M = 0.1 
model, but with the graph extended all the way back to t = 0. 

the behaviour of the unperturbed 5-tetrahedra model. 
We begin by examining the expansion rate of the uni¬ 
verse’s volume dU/dt against the volume U itself, a re¬ 
lation which has been plotted in Figure 0 The volume 
of the unperturbed universe is given by ED and its ex¬ 
pansion rate by (1321) . To first order in the perturbative 
expansion, the volume of a perturbed universe is given 

by 

C/ = J_ Z 3 + _i_, 2 (jiflW+sflW), (119) 
and the expansion rate by 


2 sin 8 


1 + 2 cos 9 
1 + cos 9 


VI. DISCUSSION OF THE MODELS 

We shall now examine the behaviour of the Regge 
model just obtained, comparing the behaviour for var¬ 
ious mass perturbations against each other and against 


dU 
tit 


■.VI 


. , -^Z'(2<H(°> + 3<H (1) ) 
+ ^j=l 2 (2 6i m +35P^ , 


( 120 ) 
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FIG. 7. Plots of Si (0) (t) (top) and 51^(t) (bottom) against t. 
The graphs have been extended all the way back to t = 0. 


FIG. 8. Plots of 5l^°\t) (top) and Sl^(t) (bottom) against t. 
The graphs have been extended all the way back to t = 0. 


TABLE II. The numerical values for 1, SI 1 ' 0 ' 1 , and 51^ at t = 

Tnax- 


SM 

M 

i 

o' 

si w 

0.001 

5.16191 x 10 -8 

-1.5465 x 10” 11 

0 

0.01 

5.16191 x 10 -8 

-1.5465 x 10 -10 

0 

0.1 

5.16191 x 10 -8 

-1.5465 x 10 -9 

0 


where 81 would be given by the explicit time-derivative 
of 81^. Across all perturbations in mass, the evolution 
of the universe’s volume is very stable, indeed closely re¬ 
sembling the evolution of the unperturbed universe. The 
effect of increasing the perturbation 8M/M is for the uni¬ 
verse to attain larger volumes and faster expansion rates. 
All graphs have been truncated on the left at the moment 
the struts turn null. However, we note that these models 
actually remain well-behaved past this point all the way 
back to t = 0, as Figure [6] shows for the SM/M = 0.1 
model. 

Figure [7] shows the behaviour of Sl^ and 51 W as func¬ 
tions of time; these graphs have also been extended all 
the way back to t = 0. They reveal that 81^ and SI W are 
indeed well-behaved near t = T max and approach zero as 
t — > T max , which is required for a moment of time symme¬ 
try. They also reveal that 51^ and 81 ^ start from zero 
as well at t = 0. Finally, they show that the absolute 
magnitudes of the graphs increase with 8M/M. 

We have computed the numerical values for l , 5l^°\ 


and SI at t = T max to verify that they are indeed zero; 
the results are given in Table [TT] As far as the computer 
is concerned, 51 ^ is exactly zero for all 5M/M. All other 
quantities are infinitesimally small, and the difference 
from zero can be attributed to numerical error. Thus, 
we can conclude that l , and SI W are indeed all zero 
at t = Tmax, and hence t = T max is a moment of time 
symmetry, as required for our choice of initial conditions 
to be valid. 

Figure [5] shows the behaviour of the length perturba¬ 
tions SI and 81 W as functions of time; these graphs 
have also been extended back to t = 0. We see that the 
magnitudes of these graphs also increase with 5M/M, 
which is consistent with the models’ attaining larger vol¬ 
umes in Figure [5] as the perturbation is increased. We 
also see that 51^ and 51 ^ are always well-behaved and 
non-zero. The latter fact implies that there is never a 
moment when the tetrahedra have edges of zero length. 
Rather, there is actually a moment when the tetrahedra 
are equilateral; this happens at t = 0, because 9 is then 
zero, and from (I118D . it therefore follows 81^ and SI W 
are equal. This suggests that our choice of a = /? = | 
for co-moving particle trajectories at the end of Section 
HVBI was appropriate. Finally, we see that the length per¬ 
turbations decrease with time, consistent with the fact 
that both 81^ and 51^ are always negative, as Figure [T] 
shows. 

Before closing, we wish to comment on our choice of 
approach to the perturbed Regge models. These models 
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were based on a specific triangulation of the CW skele¬ 
ton where all diagonals on Type II quadrilateral hinges 
terminated at E'. This triangulation had the virtue of 
simplicity, but there were other equally valid but in¬ 
equivalent triangulations we could have worked with. In 
terms of the triangulation algorithm described in Sec¬ 
tion IIV Al we used a labelling of the Cauchy surface’s 
five vertices such that the common apex to the Type II 
tetrahedra was ordered last. By shifting the order of that 
vertex’s label, one could generate an alternative triangu¬ 
lation that would certainly not be equivalent to the one 
we used. For instance, if the orders of vertices A and E 
were swapped, then the diagonals on all hinges attached 
to strut AA' would now terminate on A! instead; none 
would terminate on A. Thus, compared to the original 
triangulation, some of the diagonals on Type II quadri¬ 
lateral hinges would get swapped, and the new pair of 
triangular hinges that result would be geometrically dif¬ 
ferent from the original pair. It remains to be seen how 
many distinct triangulations there are, although we know 
there must be at least five, as there are five possible or¬ 
derings of the label for the Type II tetrahedra’s common 
apex; given a specific ordering for this vertex, it remains 
to be seen whether permutations of the other vertices’ 
ordering would generate inequivalent triangulations or 
not. More importantly, it remains to be seen whether 
these alternative triangulations would lead to the same 
Regge equations or to something new. For example, it 
is possible that in the continuum time limit, the alterna¬ 
tive triangulations would still reduce to the same Regge 
equations at leading order in dt. Even if we had a new 
set of Regge equations, the solutions of both this new 
set and the original set would be equally valid as both 
of them would satisfy the global Regge equation, as dis¬ 
cussed in Section HV Al We shall leave a more thorough 
investigation of these alternative models to future study. 
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APPENDIX A: REGULAR LATTICES IN 
3-SPACES OF CONSTANT CURVATURE 

In this appendix, we shall list all possible lattices that 
cover 3-spaces of constant curvature with a single regu¬ 
lar polyhedral cell. The cell is tiled to completely cover 
the 3-space without any gaps or overlaps. This tes¬ 
sellation problem has been thoroughly studied by Cox- 
eter [221. Clifton and Ferreira 0 have succinctly sum¬ 
marised Coxeter’s results relevant to our discussion, and 


TABLE III. All possible lattices obtained by tessellating 3- 
spaces of constant curvature with a single regular polyhedron. 
We note that the second column, which indicates how many 
cells meet at any lattice edge, effectively determines the lat¬ 
tice’s structure. 


Elementary 
cell shape 

Number of 
cells at a 
lattice edge 

Background 

curvature 

Total cells 
in lattice 

tetrahedron 

3 

+ 

5 

cube 

3 

+ 

8 

tetrahedron 

4 

+ 

16 

octahedron 

3 

+ 

24 

dodecahedron 

3 

+ 

120 

tetrahedron 

5 

+ 

600 

cube 

4 

0 

oo 

cube 

5 

- 

oo 

dodecahedron 

4 

- 

(X) 

dodecahedron 

5 

- 

oo 

icosahedron 

3 

- 

oo 


we have presented their summary in Table uni 

The Coxeter lattices are the only possible lattices that 
use a single regular polytope as its elementary cell. How¬ 
ever, if we allow for elementary cells that are not regular 
polytopes, then further regular lattices are possible. For 
instance, one can obtain a new lattice from the closed 
600-tetrahedra lattice by using the centres of the origi¬ 
nal lattice’s triangles as the new cell centres; one would 
then partition out new cells by erecting new boundaries 
between pairs of nearest-neighbouring triangular centres. 
Each tetrahedra has four triangles, but each triangle is 
shared between two tetrahedra, so the 600-tetrahedra lat¬ 
tice has a total of 1200 triangles; thus this new lattice has 
a total of 1200 cells. Clearly, this cannot be a Coxeter 
lattice, and its cell therefore cannot be a regular poly¬ 
tope; yet this new lattice is regular because the triangles 
are distributed in a regular manner. For an analogous sit¬ 
uation, consider a lattice tessellating flat 2-dimensional 
space with equilateral triangles, as shown in Figure 0 
and a new lattice with cells centred on the mid-points of 
the original edges. This new lattice is still regular, but 
its cells do not correspond to regular polytopes: their 
internal angles are not equal, even though their edges 
are; thus the cells are actually rhombi. Since the internal 
angles are not equal, the vertices themselves are not iden¬ 
tical either, with three new cells meeting at vertices like 
the one marked A and six at vertices like the one marked 
B. Therefore, such cells would not have as high a rota¬ 
tional symmetry as its regular counterpart, the square; 
the rhombus has only an order 2 rotational symmetry, 
while the square has order 4. For such reasons, the Cox¬ 
eter lattices are the lattices with the highest symmetry. 

For the closed lattices, if one constructs new lattices 
using cells centred on the original lattice vertices, it ap¬ 
pears these new lattices form Coxeter lattices as well. We 
shall refer to these new lattices as dual lattices. The 5- 
tetrahedra lattice is dual to itself, as is the 24-octahedra 
lattice. The 8-cube and 16-tetrahedra lattices are du- 
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FIG. 9. A lattice of equilateral triangles tessellating flat 2- 
dimensional space with a new lattice, derived from the orig¬ 
inal, centred on the mid-points of the triangular edges; the 
triangular lattice has been drawn in solid lines and the new 
lattice in dashed lines. The new lattice cells do not corre¬ 
spond to regular polytopes but are instead rhombi; although 
their edge-lengths are equal, their interior angles are not. For 
the new cell marked ABCD, the angles at the centres of trian¬ 
gles, marked A and C, are 120°, while the angles at vertices 
of triangles, marked B and D, are 60°. These angles imply 
that three dual cells would meet at A and C, and six at B 
and D; thus the cell vertices are not identical either. 


als of each other, as are the 120-dodecahedra and 600- 
tetrahedra lattices. However, new lattices using cells cen¬ 
tred on the mid-points of edges or the centres of faces are 
not Coxeter lattices in general. 


APPENDIX B: VARIATION OF PARTICLE 
PATH-LENGTHS WITH RESPECT TO THE 
STRUTS 


the struts, we can greatly simplify our co-ordinate sys¬ 
tem if we first impose the symmetries on all other edges, 
constraining them to be 4°4 4+u l^\ 4+u df E , df E f, 
df B , dfff accordingly; this is permissible because when 
we locally vary with respect to one edge, all other edges 
must be held constant. We then calculate Si in this new 
co-ordinate system, differentiate it with respect to each 
of mf, mf , mf, mf, and then impose the relevant strut- 
length constraints on mf, mf, mf, mf. 

Once again, the 4-block has been triangulated in the 
manner described in Sections ITT1 and II V Ai and this intro¬ 
duces the diagonals AE', BE', CE' , which have length 
df E , and the diagonals AB', AC', and BC', which have 
length df B '. 

We shall now construct our new co-ordinate system. 
We can freely fix the co-ordinates of the upper tetrahe¬ 
dron’s vertices to be 


A' = 

B' = 

C' = 


7(1) 7(1) 

n +1 L i +1 


2 v / 3 ’ 


0, 0 


l (1) 

L i +1 

2 ' 


Z (1) 

L i +1 

2 v / 3 


0,0 , 


;(i) 


0 , 0 


’ V3’ 


E' — (0,0, hi+i, 0). 


(Bl) 


Since the upper tetrahedron’s co-ordinates are fixed, the 
dependence on the strut-lengths must appear in the lower 
tetrahedron’s co-ordinates. 

Vertex A is constrained by the lengths 


AA! 


= mf. 


AB' = AC' =d: 


,AB' 


AE' =d 


AE' 


As A is equidistant to B' and C', its co-ordinates will 
have the form 


A = 



CA ’ 



In this appendix, we shall explain the derivation of 
results 4ZZD, (GB, and (ED, which give the local vari¬ 
ation of the particle path-lengths with respect to each 
strut-length. Let us consider the length Si of an arbi¬ 
trary particle through a 4-block; we shall work with the 
4-block of a Type II tetrahedron because it is more gen¬ 
eral; it can easily be reduced to the Type I case by setting 
= 0, if' 1 = lf\ and 4+\ = 4+\- In order to vary s* 
with respect to each strut locally, we need to express it 
first in terms of the lengths of all four struts. We de¬ 
note the four strut-lengths by mf, mf, mf, and mf, 
with the superscript labelling the lower vertex to which 
the strut is attached. Our approach will be to use a 
new co-ordinate system for the 4-block such that the co¬ 
ordinates of the vertices are given in terms of the lengths 
of the edges, including the tetrahedral edges, the diago¬ 
nals, and the struts. As we are only interested in varying 


Using our new co-ordinates for the vertices, we can cal¬ 
culate the edge-lengths above in terms of a a, ca, and dA- 
This leads to the equations 

) 2 = \ (4+i ~ oa) + ^ (4+i - oa) +c\- d\, 

{ d t B ') 2 = \ (4?i + «a) + ^ (44\ - a A )\c\ ~ d 2 A, 

(df E ’) = + (h i+1 - c A f- d\. 

Since we are only interested in the first derivative of Si 
with respect to the strut-lengths, we need only determine 
a A, ca, and dA and similar quantities to first order in 
Smf, Smf, Smf, and 5mf. So by expressing mf as 
mf ~ mf A + Smf and making use of (l53l) to (1551) , we 
can solve the above system of equations to first order in 
Smf, obtaining 
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r^AA' 

a A ~ Z, 1 — 2 V <5to. 


Z (1) 

‘i+i 


A 
i i 


61 


(i) 


1 m 


AA' 


CA 


d A 


r cosh ipi + sinhV’j Sti + 1 — Smf, 

2y6 3 fii+i 


,(D 


(5Z. 


(i) 


2-\/6 


sinh ipi + cosh i/ji Sti + 


1 m AA' ^6 cosh ^ + sinh ^ St * “ 


3 hj 


»+i 


z u) 

2 F rr + 1 

i i+i 


«i+i 


id) 


sinh ipi + cosh ^ 


■ 5m{ 


Vertex £? is constrained by the lengths 


BB' 


BC' 


= d 


AB' 


BE' 


= d 


AE' 




= 1 


(i) 


and we shall express its co-ordinates in the form 


B = ( -a B , --^=b B , cb, ~id B 


Furthermore, based on symmetries and the co-ordinates 
just obtained for A, we know what the co-ordinates of B 
should be when mf = mf A , so we can express a B , b B , 
cb, and d B as 

a B ~ 1 7 + 8a B , 
bs^l^+Sbs, 
si W 

cb ~ cosh?/)j + sinh i)>i Sti + 6c B , 

Si w 

d B « — l —pz sinhV’i + cosh tbi Sti + 8d B , 

2^6 

where Sa B , Sb B , Sc B , 6d B are linear in Smf and 8mf. 
We can use the new co-ordinates to express the lengths 
above in terms of as, b B , c B , and d B , yielding 


i2 

B ) 


( m ? ) 2 = J (4+1 ^ as) + ^ (*i+i - + c| - d 

( df B ) = -a B + - ^ 4+1 + - b B ^j + c 2 B - d 2 B , 

( d t E ) = \ a B + -^b 2 B + (h i+1 -c B ) 2 -d? B , 

=^{a A + a B f + J 2 ( & b - “a) 2 
+ (cb - ca ) 2 - (d B ~ dA ) 2 ■ 

Then by taking mf ss mf A + Smf and making use of 
(1531) to (1561) , we match the first order terms and solve to 


I- 

obtain 

Sa B sa 
Sb B w 
Sc B 


rnAA 

2^_ Smf, 




‘i+l 


- 2 


m. 


AA' 


Z (1) 

b+i 


(<5mf +2 Smf) 


1 m. 


AA’ 


3 hi -)_i 


Sd B « - -mi 


AA' 


(<5mp + 2 Smf ) 
Sl {1) 


2\/& 


sinh il>i + cosh tfi 


1- 


cosh i/ji +sinh V’i Sti 


h 


■i+l 


(Smf+ 2 Smf) 


/C 1 ) 


- (Smf - dm?) 


Z (1) 

b+i 


Vertex C is constrained by the lengths 


(B2) 


CC' 


jd = Ed 


=e\ 


CE' 


= d. 


AE ' 


and we shall express its co-ordinates in the form 


C = ^ac, ^&c, Cc ) ~ L dcJ ■ 

We also know what the co-ordinates of C should be when 
to? = mf A , so we can express ac , be , cc, and dc as 

ac ~dac, 

b c » Zf } + Sb c , 

sf 1 ^ 

Cc ~ cosh'i/ij + sinh^ Sti + Sec, 

2V6 


dc 


<5Z. 


(i) 


2y/6 


sinh'i/ij + cosh^>j Sti + Sdc, 


where <5ac, Sbc, Sec, Sdc are linear in Smf, Smf, and 
Smf. Following a similar method to that of the previous 
two vertices, we find that 


Sac 


mf A ' 

7^ 

b+i 


(Smf - Smf) , 
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Sb c 

Sec 

Sd c 


^5 

7« 


(Smf + Smf) 


*i +1 


1 to; 


3 h. 


AA' 


i +1 


-to A j4 ' 
3 1 


(6mf + <5mf + 3 Smf ) 




(i) 


2^6 


sinh ipi + cosh ipi Sti 


-l 


<5i 


(i) 


cosh 0j + sinh ipi Sti 


- 1 


^i+i 


x (Smf + Smf + 3 Smf ) 


/ (1) 

+ -fry + Smf) 

i+1 


Finally, vertex E is constrained by the lengths 


EE' 


Ajf = irk = ct 


= i 


(o) 


J 


and we shall express its co-ordinates in the form 
E = ( d E , bE, Ce, —uIe) ■ 


We also know what the co-ordinates of E should be when 
mf = mf E , so we can express a E , b e, ce, and d E as 


a E 

bE 


ce 


d E 


Sa E , 
Sb E , 


s_i f 

2 v / 6 

61^ 


2^6 


cosh ^ + sinh0j Sti + hi coshc^ + Sc E , 
sinh0j + cosh ipi Sti + hi sinheti + Sd E , 


where on is a yet to be determined boost parameter and 
Sa E , Sb E , Sc E , Sd E are linear in Smf, Smf, Smf, and 
Smf. Unlike the other struts, mf has the perturbative 
expansion mf ft mf E + Smf. We can follow a similar 
method to that of the previous vertices to solve the equa¬ 
tions above for a*, Sa E , Sb E , Sc E , Sd E - By matching the 
zeroth order terms, we deduce that 


Oti = Ipi- 


(B3) 


Next, by matching the first order terms and then solving, 
we find that 


Sa 


■E 


Sb 


] E 


TO 


AA' 


l {1) 

L i +1 


- Smf — 


1 TO 


AA' 


h- 

-fry + 2—^— cosh ipi 
/W h. 

l i +1 


1 mf A 


3 iW 

hi sinh ipi 


T'i+l 


Sm: 


B 


Q 7 ( 1 ) X/C 1 ) 

h V/e sinh ipi + cosh 0* Sti 


l W 

L i+1 


fil (!) 

cosh ipi + sinh 0* St 


(3 Smf — Smf) — 2 [ 1 — 


1 mf A 


1 + 


hi sinh ipi 


sinh ipi + cosh ipi Sti 


(Smf + Smf) - 


2 TO 


AA' 



\/3 lV K+i £L s i n h + CO sh xpi Sti 


Smf 


Sc E ~ — 


Sdf 


mf E sinh ipi 
Sti + hi+i sinh0, 

mf E cosh ipi 
SU + hi -)-i sinh0, 


■ Smf 


1 m 


AA' 


3 hi 


+1 


1 mf A 


1 + 


hi sinh ipi 
C 1 ) 

Vfj g Sinh ipi + cosh ipi SU 


(Smf + Smf + Smf) 


Smf + f^i 


3 h 


H+l 


fil (!) 

cosh ipi + sinh ipi Sti + hi cosh ipi - h i+1 
(Smf + Smf + Smf) ^-. 


—b= sinh ipi + cosh ipi Sti 


Using this new co-ordinate system, we shall now vary 
the particle’s path-length with respect to each of the 
struts. Based on the particle’s position given by <53 
and its counterpart on Cauchy surface E i+1 , the particle 
should follow a trajectory Si given by 

Si = i (JA 1 + ~BB' + CC' + EE'^j . 


We note that since s, is just a linear combination of 
the four strut vectors, it will always be time-like if all 
four strut vectors are time-like. For each vertex X = 
A, B, C, E, let us express the perturbative expansion of 

the corresponding strut vector XX' as 


XX' ft m x + Sm x H-, 


where m x denotes a vector corresponding to the zeroth 
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order component of XX' and Sm x denotes a vector cor- responding to the component of XX' that is first order in 

5m ?, 8m ?, 5m? , and 5m? . Then the trajectory length 
Si can be expressed as the perturbative expansion 


Si = s 0 


= So 


So 

l s o| 


• Ss 


d^? dm ' 


dsi 

dm? 


5m, 


dsj 

dm? 


B 1 w ° l 5m? + d S ? 5m? 


dm? 


where Sq denotes the vector 


and 5s the vector 


s 0 = - (m A + m B + m c + m E ), 


5s = - (5m A + 5m B + 5m c + 5m E ). 

In the new co-ordinate system, it can be shown that so 
has co-ordinates 


/ 1 / Si?) \ L (Sl K . 

so = I 0, 0, - (h i+ 1 — hi cosh'! p? — I coshipi + sinh^j 5t% J , — hi smhipi + t I sinh^j + cosh^i 5U 


/(!) 


Since the first two co-ordinates are zero, then to calculate t^t • 5s and hence dsi/dm? , we need only the third and 

fourth co-ordinates of 5s. The third co-ordinate is 


1 m. 


3 h. 


AA' 


\ hi sinh fa 


i+1 


/CD 


sinh ipi + cosh ipi 5t i/ 


( 5m? + 5mf + 5m?) + 


C\ , 2 m f E ' sinh ipi 


5ti + hi + \ sinh'0, 


5m?, 


and the fourth is 


(i) 


^coshfc + smhi/fiSti-hi+i + lhiCoshipi ± m ? E ' cosh ipi 

-(om* + 5mi + dm, j — ^- t-——- 8m 6 . 


3 h. 


i+l 


51 


(i) 


sinh ipi + cosh ipi 5ti 


5ti + h i+ 1 sinh ipi 


I- 

We can then obtain dsi/dm? for each X by reading off the factor multiplying 5m? in y|^j ■ 5s, and we find that 

I 


dsi dsi 
dm? dm? 


dsj 

dm? 


dsi 

dm? 


1 m? A> ( \ hi sinh ipi 

^ I S ° I \ sinh ipi + cosh ipi 5ti 

1 m? E ' \ h i+ 1 sinh ipi + 5ti 
4 | s 0 1 h i+1 sinh ipi + 5ti 


(B4) 

(B5) 





























33 


where 


ffi (siv M 

2 

1 Sl {1) 

2 “ 

- (hi+ 1 — hi cosh xjji) — 1 ~^j= cosh^ + sinh 7 /^ oti J 


— hi sinh ibi H- L —= sinh ibi + cosh ibi oti 

4 2^/6 



s o = 


Taking the continuum time limit, we have that 

|so| —> so dt + 0(dt 2 3 4 5 ) 

2 



- ( \ht}> + l 


dt - 


and in this limit, it follows that 
ds ds ds 1 ™ AA 


0(dt 2 ) , 

(B7) 


dm A dm B 


dm c 

ds 


1 TO 

4 


So 


-hip+ 1 ) + 0(dt), 


1 fn EE ' jhip + 1 


<9m £ 4 s 0 hip +1 

which are relations C3 and (1751) as required. 


+ O(dt), 


(B6) 

I- 

Finally, to obtain the equivalent results for the par¬ 
ticle in the Type I tetrahedron, we simply set ^ to be 
zero, replace vertex E with vertex D, and replace the 
lengths lf \ Z-°\, and m EE ' with and m AA ', 

respectively. Since the mass in this tetrahedron is the 
perturbed mass, the path-length is now s'. We then find 
that 


ds' 


ds' 


ds' 


ds' 


l m. 

4 


AA' 


dmf dmf dmf dmf 4 5ti 
and in the continuum time limit, this becomes 
ds' ds' ds' ds' 


(B8) 


dm A dm B dmp 

which is relation m as required. 


_= —— rh AA ' 

dm D 4 
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